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I.  INTRODUCTION 


The  Mediterranean  Sea  is  an  important  area  to  research  for  the  navy’s  tactical 
issues,  especially  given  its  historical  importance  and  its  position  surrounded  by  NATO 
countries.  Chapter  I  defines  the  basic  structure  of  the  Mediterranean  Sea  and  the  objective 
of  this  thesis  to  highlight  the  relevance  of  the  inter-annual  variability  of  sound 
propagation. 

A.  THE  MEDITERRANEAN  SEA 

The  word  Mediterranean  derives  from  the  Latin  medius  terrae,  meaning  “inland.” 
Surrounded  by  the  continents  of  Europe,  Asia,  and  Africa,  the  Mediterranean  is  an 
intercontinental  sea  extending  from  the  Atlantic  Ocean  to  the  Middle  East,  separating 
Europe  from  Africa  (Tanabe  2013).  Its  major  subdivisions  from  west  to  east  include  the 
Balearic,  Ligurian,  Tyrrhenian,  Adriatic,  Ionian,  and  Aegean  Seas.  Figure  1  shows 
countries  that  border  the  Mediterranean  Sea. 
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Figure  1.  The  Mediterranean  Sea,  Bordering  Countries,  and  Subdivisions. 

Adapted  from  WorldAtlas  (2016). 
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B.  GEOGRAPHY 


The  overall  length  of  the  Mediterranean  Sea  basin  is  3,860  kilometers,  and  the 
utmost  width  is  1,600  kilometers.  The  Strait  of  Gibraltar  connects  the  Mediterranean  Sea 
to  the  Atlantic  Ocean.  The  depth  of  the  Strait  of  Gibraltar  varies  between  300  and  900 
meters.  Deep  basins  of  the  Mediterranean  Sea  are  well  connected  to  one  another; 
however,  the  Black  Sea  and  the  Aegean  Sea  are  connected  only  by  a  very  narrow  strait 
(Office  of  Ocean  Exploration  and  Research  2001).  The  Mediterranean  Sea  consists  of 
two  main  sub-basins:  the  Western  Mediterranean  and  the  Eastern  Mediterranean.  The 
Western  Mediterranean  Basin  lies  between  the  Strait  of  Gibraltar  and  the  Strait  of  Sicily. 
As  described  by  Alhammoud  et  al.  (2005),  the  Eastern  Mediterranean  basin  lies  to  the 
east  of  the  Strait  of  Sicily  (-150  km  wide  and  -600  m  deep),  and  consists  of  important 
regional  seas:  the  Ionian,  Levantine,  Adriatic,  and  Aegean  Seas.  The  maximum  depths 
are  approximately  3400  meters  in  the  Western  Basin  and  4200  meters  in  the  Levantine 
Sea,  separated  by  the  Strait  of  Sicily,  which  is  roughly  430  meters  deep  (Reddy  2001). 

Both  basins  include  many  regional  seas  and  straits.  As  for  regional  seas,  the 
Alboran,  Balearic,  and  Ligurian  Seas  are  within  the  Western  Mediterranean.  The 
Tyrrhenian,  Ionian,  Adriatic,  and  Aegean  Seas  are  within  the  Eastern  Mediterranean. 
Figure  2  shows  individual  seas  that  make  up  the  Mediterranean. 
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Figure  2.  Individual  Seas  that  Make  Up  the  Mediterranean  Sea.  Adapted 

from  World  Facts  Inc.  (2016). 


C.  THESIS  OBJECTIVE 

The  Mediterranean  Sea  has  been  extensively  studied  as  this  region  is  of 
operational  significance  to  the  U.S.  Navy.  Effective  undersea  warfare  (USW)  operations 
such  as  anti-submarine  warfare  (ASW),  undersea  warfare  (USW),  and  mine 
countermeasures  (MCM)  require  a  clear  understanding  of  how  sound  propagates  through 
the  water.  Therefore,  the  North  Atlantic  Treaty  Organization  (NATO)  gives  priority  to 
improved  understanding  of  ocean  acoustics  as  a  means  to  maintain  a  tactical  advantage  in 
the  Mediterranean  Sea. 

In  recent  years,  NATO-Russian  tensions  have  risen,  particularly  in  the  Eastern 
Mediterranean,  due  to  the  development  of  a  Russian  naval  base  in  Syria’s  Tartus  port. 
Because  Turkey  has  the  longest  coastline  in  the  Eastern  Mediterranean  Sea  among  NATO 
members,  those  tensions  put  Turkey  in  a  more  delicate  situation. 

For  this  reason,  depth  characterization  of  acoustic  transmission  loss  in  the 
Mediterranean  Sea  is  of  vital  importance  and  requires  more  in-depth  study.  This  thesis 
aims  to  acquire  operationally  important  results  for  ASW  by  performing  a  detailed 

analysis  of  variations  in  temperature,  salinity,  and  sound  speed  in  the  region  to  examine 
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inter-annual  variability  of  sound  propagation.  This  study  focuses  on  data  collected  in  the 
Mediterranean  Sea  from  January  1960  to  December  2014. 


D.  THESIS  ORGANIZATION 

Some  significant  aspects  of  the  Mediterranean  Sea  oceanography  are  explained  in 
Chapter  II  including  external  forcing,  thermohaline  circulation,  and  inter-annual 
variability.  The  data,  which  include  bathymetry  and  sediment  databases,  geoacoustics 
parameters,  the  synoptic  monthly  gridded-world  ocean  database  (SMG-WOD),  and  the 
generalized  digital  environmental  model  (GDEM),  are  described  in  Chapter  III.  Chapter 
IV  contains  the  methods  used  in  this  study.  Chapter  V  describes  and  compares  seasonal 
variability  of  spice  and  sound  speed  fields  in  the  synoptic  monthly  data  with  GDEM. 
Chapter  VI  contains  discussion  of  the  interannual  variability  of  temperature,  salinity 
and  sound  speed,  and  associated  acoustic  propagation  patterns.  Chapter  VII 
summarized  research  findings  as  well  as  identifies  new  areas  for  research  within  the 
Mediterranean  Sea. 
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II.  OCEANOGRAPHY  OF  THE  MEDITERRANEAN  SEA 


The  oceanography  of  the  Mediterranean  Sea  is  unique.  It  can  be  characterized  by 
very  high  salinity,  temperature,  and  spiciness.  This  chapter  explains  the  external  forcing, 
and  its  effects  on  the  thermohaline  circulation,  water  masses,  and  sound  propagation,  and 
especially  the  intra-seasonal  and  inter-annual  variability. 

A.  EXTERNAL  FORCING 

The  Mediterranean  Sea  is  a  useful  example  of  influence  of  air-sea  interaction  or 
external  forcing  on  the  local  circulation  patterns  and  thermal  structure  over  the  water 
column.  It  is  an  evaporative  basin  with  an  estimated  freshwater  (evaporation- 
precipitation-river  runoff  =  net  evaporation)  loss  of  50  to  100  centimeter  (cm)  per  year 
(Bethoux  1979;  Bryden  and  Kinder  1991).  For  this  reason,  the  Mediterranean  Sea  is  a 
concentration  basin  on  average. 

The  Mediterranean  Sea  region  has  very  hot  and  dry  summers,  which  results  in  a 
high  rate  of  water  evaporation.  During  most  of  the  year,  the  Mediterranean  Sea  is  also 
affected  by  winds  from  the  northwest  carrying  warm  and  dry  air.  Evaporation  from  these 
winds  also  results  in  a  net  loss  of  fresh  water,  as  has  dam  placement  along  the  Nile  River. 
The  Diversion  of  Nile  River  discharge  has  caused  the  salinity  of  the  Mediterranean  Sea  to 
increase  by  about  0.04  over  a  time  scale  of  40  years  (Borghini  et  al.  2014). 

B.  THERMOHALINE  CIRCULATION 

Thermohaline  circulation  is  a  term  used  to  describe  the  movement  of  water  in  the 
ocean  in  response  to  changes  of  temperature  and  salinity.  These  two  properties  affect  the 
density  of  water.  There  is  a  high-density  difference  between  the  Atlantic  Ocean  and  the 
Mediterranean  Sea,  and  this  density  variation  creates  a  thermohaline  circulation  between 
them,  which  is  the  case  for  oceans  all  over  the  world  (Topper  and  Meijer  2015).  The 
Strait  of  Gibraltar  drives  the  circulation  because  it  connects  the  Atlantic  Ocean  with  the 
Mediterranean  Sea  (Beranger  et  al.  2005).  The  water  that  inflows  from  the  Atlantic 
Ocean  is  less  dense  in  the  upper  layer  through  the  Strait  of  Gibraltar.  The  denser  Western 
Mediterranean  water  outflows  to  the  Atlantic  Ocean  under  the  inflowing  Atlantic  water. 
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This  saltier  and  denser  water  in  the  Western  Mediterranean  Sea  also  plays  an  important 
role  in  the  thermohaline  structure  in  the  Atlantic  Ocean. 


A  strong  seasonal  thermocline  occurs  in  the  Mediterranean  Sea  due  to  significant 
evaporation  in  the  summer.  Most  of  the  Mediterranean’s  annual  evaporation  surpasses 
river  runoff  and  rainfall  (Borghini  et  al.  2014).  Nevertheless,  the  Adriatic  Sea  is  one 
exception.  It  receives  a  large  amount  of  fresh  water  from  the  Po  River. 

Seasonal  variations  of  temperature  and  salinity  also  drive  the  circulation  in  the 
Mediterranean  Sea.  The  Eastern  Mediterranean  is  exposed  to  more  evaporation  than  the 
Western  Mediterranean.  Thus,  the  Modified  Atlantic  Water  enters  the  Eastern 
Mediterranean  Sea  to  compensate  the  loss  of  water  in  the  Eastern  Mediterranean  Sea. 
Conversely,  Levantine  Intermediate  Water,  which  is  saltier  and  warmer  than  the 
Modified  Atlantic  Water,  enters  the  Western  Mediterranean  Sea  as  the  deeper  layer  in  the 
Strait  of  Gibraltar  (Artale  et  al.  2002).  Figure  3  demonstrates  the  thermohaline  circulation 
in  the  Mediterranean  Sea. 


GUI. FOP  LIONS- 


Figure  3.  The  Schematic  of  Thermohaline  Circulation  in  the 

Mediterranean  Sea.  Source:  European  Ocean  Observatory 
Network  (2000). 
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c. 


WATER  MASSES 


There  are  four  fundamental  water  masses  in  the  Mediterranean  Sea:  Modified 
Atlantic  Water,  Levantine  Intermediate  Water,  East  Mediterranean  Deep  Water,  and 
Levantine  Surface  Water.  Figure  4  depicts  the  vertical  distribution  of  the  Mediterranean’s 
water  masses. 


10’W  0"  IQ’E  are  30’E  Longitude 


Figure  4.  Mediterranean  Sea  Vertical  Distribution  of  Water  Masses. 

Source:  GRID-Arendal  (2013). 
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1.  Modified  Atlantic  Water  (MAW) 

Robinson  et  al.  (2001)  show  that  the  flow  between  the  Atlantic  Ocean  and  the 
Mediterranean  Sea  depends  primarily  on  the  density  difference  between  two  important 
water  masses:  the  Modified  Atlantic  Water  (MAW)  and  the  Levantine  Intermediate 
Water  (LIW).  From  the  data,  we  found  that  the  MAW  is  less  dense  than  the  LIW.  Thus, 
the  less  dense  MAW  flows  in  the  upper  part  of  the  Strait  of  Gibraltar.  The  denser  LIW 
flows  out  in  the  deep  layer  of  the  Strait  of  Sicily  and  reaches  the  Strait  of  Gibraltar.  It 
also  flows  out  in  the  deep  layer  of  the  Strait  of  Gibraltar.  These  two  important  water 
masses  adjust  the  mass  transportation  (Robinson  et  al.  2001).  Figure  5  demonstrates 
currents  in  the  Strait  of  Gibraltar. 
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Figure  5.  Currents  in  the  Strait  of  Gibraltar. 

Source:  National  Oceanography  Centre  Southampton  (2015). 


The  Modified  Atlantic  Water  is  generally  distinguished  as  having  a  minimum 
subsurface  salinity  of  -38.6  in  depths  of  30  to  200  m.  Its  salinity  increases  and  its  depth 
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range  decreases  as  it  flows  eastward  as  a  result  of  high  evaporation  (Alhammoud  et  al. 
2005).  Figure  6  presents  the  mean  salinity  on  the  surface  from  1945  to  2014. 
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Figure  6.  Mean  Salinity  at  the  Surface  from  GDEM  database 


Figure  6  shows  that  the  Modified  Atlantic  Water  is  concentrated  in  the  south 
region  of  the  Western  Mediterranean  Sea  in  a  narrow,  long  line.  It  flows  along  the 
Algerian  Coast  to  the  Strait  of  Sicily.  Its  mean  salinity  is  about  36  psu  at  the  Strait  of 
Gibraltar  and  increases  as  it  flows  eastward. 

2.  Levantine  Intermediate  Water  (LIW) 

Another  essential  water  mass  in  the  Mediterranean  Sea  is  the  Levantine 
Intermediate  Water  (LIW).  Citing  Tchemia  1980,  Ozsoy  et  al.  1989,  and  Ozsoy  et  al. 
1993,  Alhammoud  et  al.  (2005)  provide  a  comprehensive  description  of  the  LIW: 
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“A  homogenous  water  mass  occupies  the  intermediate  depth  between  300  and  700  m, 
below  the  AW.  It  is  well  known  that  this  water  mass  is  formed  in  winter  in  the  LB  by  the 
transformation  of  AW.  The  LIW  signature  is  characterized  by  the  highest  sub-surface 
salinity  (38.8-39.0  psu)  in  the  LB.  As  LIW  flows  westward,  its  salinity  decreases  due  to 
the  mixing  with  the  adjacent  water  masses”  (p.  9). 

Figure  7  shows  the  LIW  in  August,  and  Figure  8  shows  the  LIW  in  January.  The 
density  of  the  LIW  varies  between  28.95  and  29.10 

m. 
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Figure  7. 


LIW  T-S  Diagram  in  August  from  the  GDEM  database 
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Figure  8.  LIW  T-S  Diagram  in  January  from  the  GDEM  database 


3.  Eastern  Mediterranean  Deep  Water  (EMDW) 

The  third  water  mass  affecting  the  Mediterranean,  the  Eastern  Mediterranean 
Deep  Water  (EMDW),  is  identified  as  having  salinity  of  38.7  psu  and  temperature  of 
13.8°C  at  depths  greater  than  800  m.  It  is  formed  in  the  Adriatic  and  Aegean  Seas  in  the 
Levantine  Basin  due  to  strong  cooling  (Alhammoud  et  al.  2005).  The  density  remains 

/c  0 

between  29.05  and  29.10  — .  Figure  9  shows  the  EMDW  in  August,  and  Figure  10  shows 
the  EMDW  in  January. 
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Figure  10.  Eastern  Mediterranean  Deep  Water  T-S  Diagram  in  January 
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Eastern  Mediterranean  Deep  Water  T-S  Diagram  in  August 
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4. 


Levantine  Surface  Water  (LSW) 


Another  water  mass,  the  Levantine  Surface  Water  (LSW),  is  thin,  warm,  and 
salty.  It  resides  over  the  surface  of  the  Levantine  basin.  It  forms  particularly  in  summer 
because  of  the  high  rate  of  evaporation  (Alhammoud  et  al.  2005).  It  floats  over  the 
MAW.  We  found  that  the  LSW  is  the  saltiest  and  warmest  water  mass  among  all  others 
in  the  Mediterranean  Sea  (see  Figure  11). 
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Figure  1 1 .  Levantine  Surface  Water  in  August 
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III.  DATA 


Chapter  II  explained  the  different  forces  acting  on  the  Mediterranean  Sea.  This 
chapter  shows  how  the  four  environmental  datasets  were  utilized  into  an  acoustic  model 
to  study  the  inter-annual  acoustic  variability  (Allen  and  U.S.  Navy  2012,  Chu  and  Fan 
2016b,  NAVO  2006,  NAVO  2009,  and  Porter  2011):  sediment,  bathymetry,  SMG-WOD 
and  GDEM  salinity  and  temperature  data. 

A.  SEDIMENT  DATABASE 

The  Naval  Oceanographic  Office  (NAVO)  has  four  open  sediment  databases: 
enhanced,  standard,  reduced,  and  High-Frequency  Environmental  Acoustics  (HFEVA). 
The  enhanced  database  is  the  actual  database,  which  NAVO  maintains.  The  other  three 
sediment  databases  are  rearranged  from  this  data  source. 

The  HFEVA  sediment  database  has  23  standard  types  based  on  sediment  mixture 
and  grain-size.  There  are  two  uncommon  numbers,  888  alluding  to  “no  data”  and  “999” 
alluding  to  land.  Table  1  shows  the  HFEVA  sediment  types.  See  Appendix  A  for  a 
complete  sediment  list. 

The  HFEVA  dataset  is  derived  for  high-frequency  acoustic  performance 
prediction  (NAVO  2003).  The  23  sediment  types  are  transformed  into  a  set  of  three 
geoacoustic  parameters  for  the  BELLHOP  acoustic  model. 
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Table  1.  HFEVA  Sediment  Categories.  Source:  NAVO  (2006). 


HFEVA  Standard  Sediment  Type 

HFEVA  Category 

Rough  Rock 

1 

Rock 

2 

Cobble  or  Gravel  or  Pebble 

3 

Sandy  Gravel 

4 

Very  Coarse  Sand 

5 

Muddy  Sandy 

6 

Gravel  Coarse  Sand  or  Gravelly  Sand 

7 

Gravelly  Muddy  Sand 

8 

Medium  Sand  or  Sand 

9 

Muddy  Gravel 

10 

Fine  Sand  or  Silty  Sand 

11 

Muddy  Sand 

12 

Very  Fine  Sand 

13 

Clayey  Sand 

14 

Coarse  Silt 

15 

Gravelly  Mud  or  Sandy  Silt 

16 

Medium  Silt  or  Sand-Silt-Clay 

17 

Sandy  Mud  or  Silt 

18 

Fine  Silt  or  Clayey  Silt 

19 

Sandy  Clay 

20 

Very  Fine  Silt 

21 

Silty  Clay 

22 

Clay 

23 

No  data 

888 

Land 

999 
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B.  BATHYMETRY  DATABASE 


We  extracted  bathymetry  data  for  the  Mediterranean  Sea  from  the  NAVO  Digital 
Bathymetric  Database- Variable  resolution  (DBDB-V)  version  5.4.  There  are  three 
accessible  database  levels  of  DBDB-V.  We  used  Level  0,  which  is  unclassified,  for  this 
research.  The  DBDB-V  is  a  fine-scale  database  due  to  four  resolutions  of  bathymetry 
data:  in  arc  minutes  of  2’,  1’,  0.5’,  or  0.1’.  The  two-minute  arc  grid  had  been  generated 
by  aggregating  data  from  different  freely  accessible  sources  in  2004,  and  it  is  the  main 
resolution  that  presents  complete  global  coverage  (McDonald  2016).  This  search  uses  a 
2’  resolution. 

C.  GEOACOUSTIC  PARAMETERS 

The  geoacoustic  parameters  are  obtained  to  compute  attenuation  and  sound  speed 
in  the  bottom  sediments.  There  are  three  main  geoacoustic  parameters  used  as  a  part  of 
the  present  research:  sediment  density,  the  attenuation  coefficient,  and  the  compressional 
sound  speed.  The  attenuation  coefficient  values  for  all  bottom  sediments  in  the  region  of 
interest  were  obtained  from  (Hamilton  1980).  Table  2  provides  the  other  parameters. 
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Table  2.  Geoacoustic  Parameters  Index.  Source:  Cintron  (2001). 


Bottom 

Sediment 

Composition 

Bulk 

Grain 

Size 

Index 

Long  (32  Char) 
Name 

Density 

gm/cm^ 

Sound 

Speed 

Ratio 

Wave 

Number 

Ratio 

BOULDER 

-9 

Rough  Rock 

2.5 

2.5 

0.0137 

ROCK 

-7 

Rock 

2.5 

2.5 

0.0137 

GRAVEL 

-3 

Grave,  Cobble  or 

Pebble 

2.5 

1.8 

0.0137 

-1 

Sandy  Gravel 

2.492 

1.337 

0.01705 

-0.5 

Very  Coarse  Sand 

2.401 

1.3067 

0.01667 

0.0 

Muddy  Sand  Gravel 

2.314 

1.2778 

0.01630 

0.5 

Coarse  Sand 

2.231 

1.2503 

0.01638 

1.0 

Gravelly  Muddy 

Sand 

2.151 

1.2241 

0.01645 

SAND 

1.5 

Sand  or  Medium 

Sand 

1.845 

1.1782 

0.01624 

2.0 

Muddy  Gravel 

1.615 

1.1396 

0.01610 

2.5 

Silty  Sand  or  Fine 

Sand 

1.451 

1.1073 

0.01602 

3.0 

Muddy  Sand 

1.339 

1.0800 

0.01728 

3.5 

Very  Fine  Sand 

1.268 

1.0568 

0.01875 

4.0 

Clayey  Sand 

1.224 

1.0364 

0.02019 

4.5 

Coarse  Sand 

1.195 

1.0179 

0.02158 

5.0 

Sandy  Silt 

1.169 

0.9999 

0.01261 

5.5 

Medium  Silt 

1.149 

0.9885 

0.00676 

SILT 

6.0 

Silt 

1.148 

0.9873 

0.00386 

6.5 

Fine  Silt 

1.147 

0.9861 

0.00306 

MUD 

7.0 

Sandy  Clay 

1.147 

0.9849 

0.00242 

7.5 

Very  Fine  Silt 

1.146 

0.9837 

0.00194 

8.0 

Silty  Clay 

1.145 

0.9824 

0.00163 

CLAY 

9.0 

Clay 

1.145 

0.9800 

0.00148 

10.0 

1.145 

0.9800 

0.00148 
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D.  GDEM 


The  GDEM  provides  monthly  climatological  fields  of  temperature  and  salinity. 
Carnes  (2004)  states,  “Development  of  GDEM  at  the  Naval  Oceanographic  Office  began 
in  1975  and  culminated  in  the  first  release  to  the  Navy  community  in  1984.  The  first 
release  contained  only  the  North  Atlantic  Region,  but  by  1991  most  of  the  world’s  oceans 
were  included.”  The  early  form  of  the  GDEM  had  a  few  restrictions  in  regions  where  the 
depth  was  shallower  than  100  meters.  In  this  manner,  the  main  adaptation  of  GDEM  did 
not  make  a  difference  between  the  shallow  and  waterfront  waters.  At  that  point,  GDEM- 
V  2.6  was  enhanced  by  another  procedure  to  empower  to  deliver  vertical  profiles  in 
regions  where  bottom  depths  are  5  meters.  Although  GDEM-V  2.6  demonstrated  better 
results  in  shallow  waters,  it  had  a  few  problems  capturing  salinity  and  temperature 
vertical  profiles  in  some  region  of  the  world’s  oceans. 

GDEM-V  3.0,  the  current  working  version,  has  been  utilized  as  a  part  of  the 
present  research.  GDEM3  was  derived  from  temperature  and  salinity  profiles  extracted 
from  the  Master  Oceanographic  Observational  Data  Set  (MOODS)  in  1995.  Staff  at  the 
Naval  Research  Laboratory  (NRL)  altered  the  profiles  and  afterward  used  them  to  build 
the  climatology  of  the  Modular  Ocean  Data  Assimilation  System  (MODAS)  by  Fox  et  al. 
(2001-2002).  Carnes  states, 

GDEM  profiles  were  interpolated  by  the  method  of  the  weighted  parabolas 
to  create  78  standard  depths  starting  from  the  surface  0m  and  going  down 
to  6600m.  It  has  a  horizontal  resolution  of  0.25.  The  gridded  of  the 
temperature,  salinity,  temperature  variance,  and  salinity  variance  was 
performed  separately  at  each  of  the  78  depth  levels. 

However,  GDEM  is  not  an  open  source  as  it  is  open  only  to  the  U.S.  military  and  its 
members. 
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E.  SMG-WOD 

Four-dimensional  (4D)  SMG-WOD  was  produced  from  the  NOAA  WOD13  at 
NPS  using  the  optimal  spectral  decomposition  (OSD)  method  Chu  and  Fan  (2016),  and 
provides  4D  world  ocean  gridded  temperature  and  salinity  data,  from  January  1945  to 
December  2014,  its  spatial  bounding  box  coordinates  are  limited  to  180E-180W  and 
80N-80S. 

The  WOD13  database  contains  a  “collection  of  scientifically  quality-controlled 
ocean  profde  and  plankton  data  that  includes  measurements  of  temperature,  salinity, 
oxygen,  phosphate,  nitrate,  silicate,  chlorophyll,  alkalinity,  pH,  pC02,  TC02,  Tritium, 
A13Carbon,  A14Carbon,  A180xygen,  Freons,  Helium,  A3Helium,  Neon,  and  plankton” 
(Boyer  et  al.  2013).  Table  3  provides  the  WOD13  datasets  and  the  instrument  sources  that 
contributed  data. 

As  a  part  of  ocean  investigation,  we  utilized  the  OSD-technique,  which  proved 
fruitful  in  drawing  out  an  ocean  phenomenon  that  had  not  documented  previously  (see 
Chu  et  al.  2005;  Chu  et  al.  2007;  and  Chu  2011)  and  might  be  utilized  as  part  of  ocean 
data  assimilation  (Chu  et  al.  2015).  This  method  makes  use  of  the  geography  and  lateral 
boundary  conditions  but  is  different  from  past  interpolation  methods,  for  instance, 
optimal  interpolation  (OI)  because  it  does  not  require  a  foundation  error  covariance 
matrix  (Chu  and  Fan  2016a).  Table  4  presents  the  vertical  levels  of  the  SMG-WOD  data 
(Chu  et  al.  2003a;  2003b;  2004;  and  2015),  which  is  now  publicly  accessible  at  the 
National  Oceanic  and  Atmospheric  Association  (NOAA)  and  the  National  Centers  for 
Environmental  Information  (NCEI)  website  (Chu  and  Fan  2016).  The  SMG-WOD  data 
for  the  Mediterranean  Sea  region  has  higher  horizontal  (0.25°)  and  vertical  (57  levels) 
resolution,  and  covers  the  period  1960-2014. 
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Table  3.  Datasets  and  the  Instruments  Sources  Contribute  Data. 

Source:  Boyer  et  al.  (2013). 


DATASET 

SOURCE 

OSD 

Bottle,  low-resolution  Conductivity-Temperature-Depth  (CTD),  low- 

resolution 

XCTD  data,  and  plankton  data 

CTD 

High-resolution  Conductivity-Temperature-Depth  (CTD)  data  and  high- 

resolution  XCDT  data 

MBT 

Mechanical  Bathythermograph  (MBT)  data,  DBT,  micro-BT 

XBT 

Expendable  (XBT)  data 

SUR 

Surface  only  data  ( bucket,  thermosalinograph) 

APB 

Autonomous  Pinniped  Bathythermograph  -  Time-  Temperature-Depth 
recorders  attached  to  elephant  seals 

MRB 

Moored  buoy  data  from  TAO  (Tropical  Atmosphere-Ocean),  PIRATA 
(moored  array  in  the  tropical  Atlantic).  MARNET,  and  TRITON  (Japan- 
JAMSTEC) 

PFL 

Profiling  float  data 

DRB 

Drifting  buoy  data  from  surface  buoys  with  thermistor  chains 

UOR 

Undulating  Oceanographic  Recorder  data  from  a 
Conductivity/Temperature/Depth  probe  mounted  on  a  towed  undulating 

vehicle 

GLD 

Glider  data 
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-WOD  Depth  Layers.  Source:  Rodriguez  (2016). 


IV.  METHODS 


Chapter  IV  explains  the  methodology  including  the  composite  analysis,  acoustic 
modeling  using  BELLHOP,  choice  of  reference  points  based  on  bottom  sediment  type, 
calculation  of  spiciness  and  quantifying  uncertainty  in  the  transmission  loss. 

A.  COMPOSITE  ANALYSIS 

Prom  the  SMG-WOD  gridded  dataset,  the  sound  speed  (c)  profdes  are  calculated 
from  the  (T,  S)  profdes.  At  a  specific  horizontal  grid,  it  is  represented  by  c(zi,Tk,tl), 

with  Zi  represents  the  vertical  levels;  i  =  1,  2,  ...,/,  where  I  (=  28)  is  the  total  number  of 
the  vertical  levels;  Tk  =  1960,  1961,  ...,  2014  is  the  time  sequence  in  years;  and  6  =  1,  2, 
...,  12  is  the  monthly  sequence  within  a  year. 

Two  Temporal  Means 

Before  investigating  the  annual  variation  of  c,  we  define  the  following  two 
temporal  averages: 

c(zi,tl)  =  ^-Y.c(zi’TkM  Ar  =  55  (1960  to  2014),  (4.1) 

Ax  V 

which  are  the  multiyear  averaged  mean  values,  and 

^(z,)  =  TvZ^Z-’r/)’  (4‘2) 

i=i 

which  is  the  total-time  mean. 

Mean  Seasonal  Variability 

The  annual  mean  values  deviated  from  the  total-time  mean,  c(z;  ) , 
c(zl,tl)  =  c(zl,tI)-c(zi),  (4.3) 
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lead  to  composite  features  of  the  mean  seasonal  variability. 

Interannual  Variability 
The  deviation  from  its  annual  mean  value 
c(Zi  ,Tk,tl)  =  c(Zi  (4.4) 

Is  called  the  anomaly  data,  which  is  re-arranged  into  a  IxP  matrix,  c(z,J p ) ,  p  =  1,2,  ..., 

P.  Here  P  =  660  is  the  total  number  of  months  (i.e.,  total  time  points).  The  empirical 
orthogonal  function  (EOF)  analysis  (Chu  et  al.,  1998;  Chu  2010)  was  conducted  on  the 
anomaly  data 

c(z,  ,tp)  =  YJ  pCj  (tp  )Ej  (Zi),  (4.5) 

j 

where  £j(z,)  is  the  jth  EOF  mode  with  the  unit  of  m/s;  and  PC  .(F  )  is  the  jth  principal 
component  without  unit.  Here,  j  =  1,  2,  ...,/;/  is  truncated  number  of  the  EOF  modes. 

B.  BELLHOP 

BELLHOP  was  used  in  this  present  research  in  order  to  predict  transmission  loss 
for  acoustic  popagation  in  the  Mediterranean  Sea.  Porter  (2011)  states,  “BELLHOP  can 
produce  a  variety  of  useful  outputs  including  transmission  loss,  eigenrays,  arrivals,  and 
received  time  series.  It  allows  for  range-dependence  in  the  top  and  bottom  boundaries 
(altimetry  and  bathymetry),  as  well  as  in  the  sound-speed  profile.”  The  sound  speed 
profiles  uploaded  into  BELLHOP  were  created  by  the  GDEM  and  SMG-WOD  for 
comparison.  BELLHOP  is  a  program  that  can  be  used  on  many  platforms.  In  this  present 
research,  BELLHOP  is  used  in  MATLAB. 

BELLHOP  creates  a  two-dimensional  ray  that  traces  plots  for  the  given  sound- 
speed  profiles  and  bottom  sediment.  Although  the  Mediterranean  Sea  is  not  shallow,  its 
mean  depth  falls  below  the  mean  for  all  global  oceans.  Additionally,  it  has  very  shallow 
waters  in  some  crucial  areas  for  NATO  operations. 
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There  are  many  reasons  why  we  preferred  BELLHOP  in  the  present  research. 
Lirst,  it  is  an  open-source  program  and  can  be  used  internationally.  Second,  modeling 
acoustic  transmissions  with  rays  is  a  common  method  for  studying  and  understanding 
how  sound  energy  propagates  through  a  given  sound  channel  (Porter  2011). 

Bottom  sediment  types  exhibit  different  parameters  in  the  calculation.  Each 
bottom  sediment  provides  a  different  bottom  reflection  coefficient.  In  order  to  analyze 
bottom  sediment  effects  on  sound  propagation,  different  bottom  sediments  were  taken 
into  account.  The  calculation  for  transmission  loss  is  shown  in  Equation  (4.6) 

I(r,  z) 

TL  =  -10  log^—^ 

=  -20  log^l  (db  re  1  m)  (4.6) 

There  are  three  different  ways  to  plot  in  BELLHOP.  They  are:  coherent, 
incoherent  and  semi-coherent.  The  user  selects  one  depending  on  the  level  of 
interference.  For  example,  if  a  user  wants  to  calculate  fine -patterned  interference  in  the 
transmission  loss,  the  coherent  mode  is  required.  Otherwise,  the  user  may  select  the 
incoherent  mode.  Semi-coherent  transmission  loss  is  a  combination  of  coherent  and 
incoherent.  It  calculates  some  interference  but  not  as  much  as  the  coherent  mode 
(McDonald  2016).  Figure  12  demonstrates  the  BELLHOP  structure. 
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The  left  side  of  the  figure  exhibits  all  the  inputs  that  can  be  uploaded  into  BELLHOP. 
The  right  side  of  figure  represents  model  outputs. 

Figure  12.  BELLHOP  Structure 


BELLHOP  transmission  losses  were  plotted  for  a  source  depth  of  40  meters  and 
various  receiver  depths.  The  results  are  explained  in  Chapter  V. 


C.  TREATMENT  OF  BOTTOM  SEDIMENTS 

The  effects  of  bottom  sediment  on  transmission  loss  depend  on  sound  frequencies. 
At  low  frequencies,  less  than  10  Hz,  bottom  sediments  parameters  are  required  a  few 
meters  above  the  sea  floor.  At  frequencies  higher  than  10  Hz,  all  sediment  thicknesses 
and  materials  underlying  the  seafloor  have  to  be  provided.  In  this  thesis  the  frequency  of 
3500  Hz  was  selected  to  calculate  transmission  loss  via  BELLHOP.  Therefore,  bottom 
sediments  are  needed  due  to  high  frequencies  penetration  under  the  seafloor. 
(Computational  Ocean  Acoustic,  2011). 
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The  Mediterranean  Sea  has  a  variety  of  bottom  sediment  at  different  depths. 
Although  the  Mediterranean  is  not  considered  shallow,  it  is  shallower  than  other  oceans 
with  a  four-kilometer  mean  depth.  Seafloor  properties  of  the  Mediterranean  Sea  were 
extracted  and  plotted  in  Figure  13. 
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Figure  13.  Bottom  Sediments  and  Bathymetry  of  the  Mediterranean  Sea 


Five  reference  points  were  selected  to  model  acoustic  propagation  in  the 
Mediterranean  Sea.  Although  the  Mediterranean  has  two  sub-basins,  the  Western  Basin 
and  Levantine  Basin,  some  authors  classify  the  water  masses  into  western,  central,  and 
Levantine.  Points,  labeled  A  through  E,  have  been  selected  for  these  water  masses  as 
shown  in  Table  5. 
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Table  5.  Locations  A  through  E  by  Bottom  Sediment  Types  and  Depth 


Location 

Latitude  (N) 

Longitude  (E) 

Depth  (m) 

Bottom 

Sediment 

Region 

A 

33.7794 

34.939 

1000 

clay 

Eastern 

B 

34.2399 

13.6275 

200 

very  fine  sand 

Central 

C 

32.1361 

27.9229 

1000 

sandy  mud 

Eastern 

D 

36 

18 

1000 

very  fine  silt 

Eastern 

E 

40.3946 

1,1739 

70 

sandy  mud 

Western 

D.  CALCULATION  OF  SPICINESS  AND  SOUND  SPEED 

Spiciness  is  a  useful  variable  for  characterizing  the  water  masses  at  a  given 
isopycnal  surface.  Flament  (2002)  states, 

It  is  simultaneously  a  tracer  conserved  by  isentropic  motions,  and  an 
indicator  of  diffusive  stability.  It  is  useful  for  the  combined  description  of 
interleaving  and  double-diffusive  processes  at  the  boundary  between  water 
masses. 

Potential  temperature  and  salinity  can  be  used  to  define  spiciness.  Positive  spiciness 
corresponds  to  warm,  salty  water.  Inversely,  negative  spiciness  corresponds  to  cold,  fresh 
water.  High  positive  spiciness  usually  results  from  the  lack  of  fresh  water  and  air-sea 
interaction  (Bindoff  and  McDougall  1994).  The  equation  from  Flament  (1986)  for 
computing  spiciness  is  as  follows: 

<9,S)  =  -  35 y  (4.7) 

Spiciness  in  this  thesis  was  calculated  using  Flament’ s  algorithm.  See  Appendix  B  for 
more  information  on  the  spiciness  calculation  function. 

The  Mediterranean  Sea  loses  buoyancy  and  exhibits  a  stratified  water  column  due 
to  considerable  evaporation,  especially  in  the  summer.  The  Levantine  Surface  Water 
(LSW)  forms  in  summer;  it  is  very  hot  and  salty  with  positive  spiciness.  Figure  14  shows 
the  spiciness  of  the  LSW  in  August. 

Sound  speed  depends  on  density  and  compressibility.  Density  is  also  related  to 
temperature  and  salinity  in  the  ocean.  Medwin’s  equation  to  express  sound  speed  is 
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provided  in  Equation  4.8.  There  are  three  variables:  temperature  (7)  in  Celsius,  salinity 
(5)  in  ppt,  and  depth  (z)  in  meters. 

c  =  1449.2  +  4.6T  -  0.055T2  +  0.00029T3  +  (1.34  -  0.0170(5  -  35)  +  0.016z  (4.8) 


33  30.5  39  39.5 

Salinity  (psu) 

Figure  14.  Spiciness  of  the  Levantine  Surface  Water  in  August 

Equation  (4.8)  indicates  that  when  temperature  and  salinity  increase  together, 
sound  speed  also  increases.  The  increase  of  salinity  and  temperature  corresponds  to 
elevated  spiciness.  The  vertical  structure  of  spiciness  affects  the  vertical  structure  of  the 
sound  speed. 

E.  TRANSMISSION  LOSS 

An  acoustic  signal  propagating  in  the  ocean  becomes  distorted  due  to  various  loss 
mechanisms.  Jensen  (2011)  says,  “The  standard  measure  in  underwater  acoustics  of  the 
change  in  signal  strength  with  range  is  transmission  loss  defined  as  the  ratio  in  decibels 
between  the  acoustics  intensity  7(r,z)  at  a  field  point  and  the  intensity  I0  at  1-m  distance 
from  the  source.” 
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Low  frequencies  propagate  further  than  high  frequencies  do.  Therefore,  low 
frequencies  can  be  more  critical  for  detection  of  submarines.  The  most  common  noise 
source  of  a  submarine  is  the  propeller.  Noise  produced  by  bubbles  collapsing  is  called 
cavitation.  At  low  speeds,  machine  noise  often  dominates  the  cavitation,  but  at 
intermediate  and  high  speeds,  especially  when  the  submarine  travels  at  periscope  depths, 
cavitation  often  dominates  all  other  noise.  Cavitation  depends  on  the  operational  depth 
and  speed.  Speeds  that  generate  cavitation  increase  with  depth  (Renilson,  2015). 

Researchers  express  the  active  sonar  equation  and  calculate  a  signal-to-noise  ratio 
(SNR)  of  possible  submarines  at  the  receiver.  Equation  (4.9)  provides  the  active  sonar 
equation. 

SNR(decibels)  =  SL  -  2 TL  +  TS-  ( NL  -  AG)  (4.9) 

TL  stands  for  transmission  loss,  TS  stands  for  target  strength,  NL  stands  for  the 
noise  level,  and  AG  stands  for  the  array  gain. 

For  this  thesis,  the  following  values  for  Equation  (4.9)  from  DOSITS.org  were 
assumed  as  reasonable. 

SL  =  220  dB;  the  typical  array  gain  is  20  dB;  TS  =  25  dB;  NL  (total)  is  73  dB;  and 
a  reasonable  SNR  is  22  dB.  The  equation  can  be  written  as  follows. 

TL  =  Q  *  (SL  -  SNR  (dB)  +  TS  -  (NL  -  AG))  (4.10) 

These  values  depend  on  the  type  of  submarine  and  ship,  the  speed,  and  other 
variables.  After  the  above  values  are  put  into  Equation  (4.10),  it  gives 

TL  =  (i)  *  (220  -  22  +  25  -  (73  -  20))  (4.11) 

TL=  85  dB 

Therefore,  values  of  85  dB  and  below  can  be  assumed  as  reasonable  TL  for 
detection  of  a  submarine.  The  transmission  loss  includes  sound-spreading  loss  and 
attenuation  loss.  We  use  the  BELLHOP  model  to  calculate  the  TL. 
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V.  SEASONAL  VARIABILITY 


This  chapter  discusses  seasonal  variability  of  the  Mediterranean  Sea  thermohaline 
structure  and  associated  acoustic  environment.  First,  we  compare  SMG-WOD  multi-year 
monthly  averaged  and  GDEM  monthly  mean  fields  of  seawater  spiciness  at  several 
horizons,  which  are  characteristic  of  different  water  masses.  Then  we  discuss  the 
corresponding  fields  of  sound  speed  following  by  discussion  of  sound  speed  profiles  in 
several  selected  locations.  The  final  subsection  contains  description  of  seasonal 
variability  of  acoustic  propagation  at  these  locations. 

A.  SPICINESS 

1.  SMG-WOD  MULTI-YEAR  MONTHLY  AVERAGE 

The  mean  spiciness  in  the  Mediterranean  Sea  showed  a  large  range  of  variability, 
both  seasonally  and  regionally.  SMG-WOD  multiyear  monthly  average  fields  of  spice  are 
shown  in  Figure  15  for  sea  surface,  and  in  Figure  16  for  the  depth  of  100  meters. 
Spiciness  was  characterized  by  strong  decrease  from  the  surface  to  the  100  m  depth, 
especially  in  August.  The  maximum  value  (9. 1  m3 /kg)  was  found  at  the  surface  of  the 
Eastern  Mediterranean  Sea,  decreasing  with  depth  to  about  100  meters  to  5. 6  m3 /kg  It 
also  decreased  in  the  Western  Mediterranean  Sea  from  6.7  m3 /kg  to  4.2  m3 /kg  at  about 
100  meters.  Spice  did  not  show  the  same  variations  from  the  surface  to  the  100  meters 
depth  in  February.  Maximum  value  (5. 6  m3 /kg)  is  found  at  the  surface  of  the  Eastern 
Mediterranean  Sea,  increasing  with  depth  up  to  5. 8  m3 /kg)  at  100  m.  It  also  increases  in 
the  Western  Mediterranean  Sea,  from  3.5  m3 /kg  to  3.8  m3 /kg  at  about  100  meters.  The 
main  reason  for  this  difference  between  spice  variations  between  summer  and  winter 
seasons  is  the  seasonal  variability  of  temperature  because  of  changing  air-ocean  heat 
fluxes.  The  surface  water  gets  warmer  during  the  summer  season  due  to  significant 
heating,  and  colder  in  the  winter  season  due  to  surface  cooling. 

The  layer  between  300  meters  and  700  meters  in  the  Eastern  Mediterranean  is 
occupied  by  the  Levantine  Intermediate  Water  (LIW).  Sometimes,  LIW  can  extend  to 
800  meters.  SMG-WOD  multi-year  monthly  averaged  spice  at  400  m  is  plotted  in  Figure 

17  to  analyze  the  Levantine  Intermediate  Water. 
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Figure  15.  SMG-WOD  Multi-year  Monthly  Averaged  Spiciness 

at  Zero  Meters 
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Figure  16.  SMG-WOD  Multi-year  Monthly  Averaged  Spiciness 

at  100  Meters 
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The  spice  decreased  significantly  at  the  depth  of  400  meters  relative  to  the 
surface.  Unlike  the  situations  at  the  surface,  the  minimum  values  were  found  in  August  in 
the  Eastern  Mediterranean  Sea.  At  400  meters  depth,  it  is  approximately  4.9  m3/kg.  The 
maximum  values  (5.1  m3 /kg)  were  seen  in  December.  There  were  high  spice  values 
around  the  Crete  passage  where  the  LIW  forms  in  winter  due  to  strong  wintertime 
northerlies.  That  finding  provided  a  good  idea  of  where  and  when  the  LIW  forms  (Millot 
2013).  High  salinity  and  temperatures,  or  positive  spiciness,  are  the  hallmarks  of  the 
LIW. 

SMG-WOD  multi-year  monthly  averaged  spice  was  plotted  at  1050  meters  to 
analyze  the  Eastern  Mediterranean  Deep  Water  (see  Figure  18).  The  maximum  value  was 
found  near  the  Crete  Passage  in  March  at  approximately  4.82  m3 /kg  The  minimum 
value  was  found  in  the  southern  part  of  the  Crete  Passage  in  June  at  approximately 
4.5 6  m3 /kg.  There  was  no  significant  difference  of  spice  as  expected  in  the  Eastern 
Mediterranean  Sea. 

2.  GDEM  MONTHLY  MEANS 

The  GDEM  monthly  spice  fields  were  analyzed  at  the  same  depths  was  analyzed 
to  compare  it  with  that  of  the  SMG-WOD  model.  Figures  19  and  20  show  GDEM 
monthly  spice  fields  at  surface  and  at  100  meters.  The  maximum  value  (9. 1  m3 /kg)  was 
found  at  the  surface  of  the  Eastern  Mediterranean  Sea,  decreasing  with  a  depth  of  about 
100  meters  to  5.  5  m3 /kg  in  August.  It  also  increased  in  the  Western  Mediterranean  Sea, 
from  6.3  m3 /kg  to  4.1  m3 /kg  at  100  m.  There  was  a  slight  difference  between  two 
models  in  spice  values  in  coastal  areas,  where  SMG-WOD  has  lower  resolution.  The 
SMG-WOD  corresponded  well  with  the  GDEM  for  deep-water  values.  Appendix  C 
provides  a  comprehensive  collection  of  spice  figures.  As  expected,  it  was  confirmed  that 
Eastern  Mediterranean  was  more  spicy  than  the  Western  Mediterranean  in  all  months. 
Additionally,  the  data  revealed  that  the  LWS,  which  forms  in  August,  has  the  highest 
spice  as  compared  to  other  water  masses  in  the  Mediterranean  Sea. 
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Figure  17.  SMG-WOD  Multi-year  Monthly  Averaged  Spiciness 

at  400  Meters 
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Figure  18.  SMG-WOD  Multi-year  Monthly  Averaged  Spiciness 

at  1050  Meters 
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Figure  19.  GDEM  Monthly  Spice  at  Zero  Meters 
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Figure  20.  GDEM  Monthly  Spice  at  100  Meters 
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B.  SOUND  SPEED 


Although  there  are  an  increasing  number  of  naval  operations  in  the  Mediterranean 
Sea,  sound  speed  profiles  in  the  Mediterranean  Sea  are  still  lacking.  One  of  the  purposes 
of  this  thesis  was  to  provide  a  background  on  sound  speed  profiles  because  the 
Mediterranean  exhibits  considerable  seasonal  and  regional  variations  due  to  variation  of 
temperature. 

1.  SMG-WOD  MULTI-YEAR  MONTHLY  AVERAGE 

SMG-WOD  multi-year  averaged  fields  of  sound  speed  at  0  m  and  100  m  are 
displayed  to  show  relation  between  spice  and  sound  speed.  Figure  21  shows  the  SMG- 
WOD  sound  speed  at  zero  meters;  the  maximum  values  for  spice  and  sound  speed  are 
summarized  in  Table  6.  In  the  Eastern  Basin,  the  sound  speed  was  higher  than  in  the 
Western  basin.  The  maximum  surface  sound  speed  (1545m/s)  was  found  in  the  east  part 
of  the  Eastern  Mediterranean  Sea  where  the  highest  spiciness  (9. 1  m3 /kg)  manifested  in 
the  August.  The  maximum  sound  speed  in  the  southern  part  of  the  central  Mediterranean 
Sea  was  approximately  1542  m/s  as  observed  close  to  the  coast  of  Libya.  For  this 
location,  spice  was  about  8.6  m3 /kg.  The  sound  speed  in  the  Western  Mediterranean 
Sea,  where  the  spice  is  6.7  m3 /kg,  was  about  1535  m/s.  As  seen  from  Table  6,  which 
summarizes  these  results,  spice  and  sound  speed  are  directly  proportional  to  each  other. 
Table  7  presents  maximum  values  of  sound  speed  and  spice  at  0  m  in  January,  which  are 
also  correlated. 


Table  6.  Sound  Speed  and  Spiciness  at  Zero  Meters  in  August 


Region 

Spice 

Sound  Speed 

Eastern  part  of  the  Eastern  Mediterranean  Sea 

9. 1  m3  /kg 

1545m/s 

Southern  part  of  central  Mediterranean  Sea 

8.6  m3 /kg 

1542  m/s. 

West  of  Sardinia  Island 

6.7  m3  /kg 

1535  m/s 
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Table  7.  Sound  Speed  and  Spiciness  at  0  meter  in  January 


Region 

Spice 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 6  m3 /kg 

1520  m/s 

Southern  part  of  central  Mediterranean  Sea 

5.4  m3 /kg 

1518  m/s 

West  of  Sardinia  Island 

3.7  m3  /kg 

1510  m/s 
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Figure  21. 


Figure  22  shows  the  SMG-WOD  sound  speed  at  100  meters;  the  maximum  values 
for  spice  and  sound  speed  for  August  and  January  are  summarized  in  Tables  8  and  9.  See 
Appendix  D  for  a  comprehensive  collection  of  the  SMG-WOD  multi-year  averaged 
monthly  fields  of  sound  speed  and  spice  for  deeper  layers. 


Table  8.  Sound  Speed  and  Spice  at  100  Meters  in  August 


Region 

Spiciness 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 6  m3  /kg 

1520  m/s 

Southern  part  of  central  Mediterranean  Sea 

4.6  m3 /kg 

1517  m/s 

West  of  Sardinia  Island 

4.2  m3 /kg 

1508  m/s 

Table  9.  Sound  Speed  and  Spiciness  at  100  Meters  in  January 


Region 

Spice 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 8  m3 /kg 

1522  m/s 

Southern  part  of  central  Mediterranean  Sea 

4.8  m3 /kg 

1515  m/s 

West  of  Sardinia  Island 

4  m3 /kg 

1508  m/s 
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Figure  22.  SMG-WOD  Monthly  Mean  Sound  Speed  at  100  Meters 
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2. 


GDEM  MONTHLY  FIELDS 


We  analyzed  the  monthly  sound  of  the  GDEM  at  same  depths  to  compare  the 
results  of  the  two  models.  The  GDEM  monthly  sound  speed  at  0  meters  and  100  meters 
are  plotted  in  Figures  23  and  24.  Other  results  are  provided  in  Appendix  E.  Tables  10-13 
summarize  maximum  values  for  spice  and  associated  sound  speed.  The  maximum  surface 
sound  speed  (1545  m/s)  was  found  in  the  eastern  part  of  the  Eastern  Mediterranean  Sea. 
This  is  in  agreement  with  the  SMG-WOD  multi-year  averaged  data. 


Table  10.  Sound  Speed  and  Spiciness  at  Zero  Meters  in  August 


Region 

Spiciness 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

9. 1  m3  /kg 

1545m/s 

Southern  part  of  central  Mediterranean  Sea 

8.6  m3 /kg 

1542  m/s. 

West  of  Sardinia  Island 

6.7  m3 /k,q 

1535  m/s 

Table  1 1 .  Sound  Speed  and  Spiciness  at  Zero  Meters  in  January 


Region 

Spiciness 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 7  m3 /kg 

1522  m/s 

Southern  part  of  central  Mediterranean  Sea 

5.3 /kg 

1519  m/s 

West  of  Sardinia  Island 

3.8  m3 /kg 

1510  m/s 

Table  12.  Sound  Speed  and  Spiciness  at  100  Meters  in  August 


Region 

Spiciness 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 5  m3/kiq 

1519  m/s 

Southern  part  of  central  Mediterranean  Sea 

4.7  m3  /kg 

1516  m/s. 

West  of  Sardinia  Island 

4.2  m3  /kg 

1507  m/s 

Table  13.  Sound  Speed  and  Spiciness  at  100  Meters  in  January 


Region 

Spiciness 

Sound  Speed 

East  Part  of  the  Eastern  Mediterranean  Sea 

5. 7  m3/kp 

1521  m/s 

Southern  part  of  central  Mediterranean  Sea 

4.9  m3  /kg 

1517  m/s 

West  of  Sardinia  Island 

4.1  /kg 

1503  m/s 
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Figure  23.  GDEM  Monthly  Sound  Speed  at  Zero  Meters 
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Figure  24.  GDEM  Monthly  Sound  Speed  at  100  Meters 
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Unless  pressure  dominates,  the  maximum  sound  speed  is  usually  found  in  regions 
where  maximum  spiciness  exists.  Therefore,  we  decided  that  the  spiciness  variable  could 
be  used  to  compare  the  water  masses  for  sound  speed  in  regions  where  pressure  does  not 
dominate.  In  deep  layer,  the  sound  speed  will  be  so  high.  Thus,  although  we  see  very  low 
spiciness  in  deep  layers,  the  sound  speed  will  be  so  high  due  to  pressure  effects.  In 
addition,  there  was  a  slight  difference  in  spiciness  values  of  the  GDEM  and  SMG-WOD; 
however,  the  SMG-WOD  did  not  perform  well  in  coastal  regions. 

C.  ACOUSTIC  PROPAGATIONS  AND  SOUND  SPEED  PROFILES  (SSP) 

Three  fundamental  layers  exist  in  the  ocean  according  to  sound  speed  (Figure  25). 
The  surface  layer  is  the  most  variable  layer  in  the  ocean.  It  can  be  affected  by  winds  and 
currents  among  other  variables.  It  is  also  called  a  well-mixed  layer.  The  main  thermocline 
layer  is  located  between  the  surface  layer  and  deep  isothermal  layers.  It  is  the  layer  where 
in  the  temperature  changes  rapidly  with  increasing  depth.  Under  the  main  thermocline, 
the  temperature  is  relatively  constant;  however,  sound  speed  starts  to  increase  due  to  the 
pressure  effect. 

In  this  research,  the  surface  layer  is  the  main  layer  of  the  interest.  It  is  also 
referred  to  as  seasonal  layer.  It  is  greatly  influenced  by  the  temperature.  During  the 
summer,  the  surface  waters  warm  more  quickly  than  the  deeper  waters.  Thus,  there  is  a 
sharp  negative  gradient  below  the  surface  water.  This  is  common  in  Mediterranean  Sea 
due  to  very  warm,  and  thin  Levantine  Surface  Water  that  forms  in  August. 

In  summer  seasons,  sound  channels  can  occur  below  the  main  thermocline.  At  the 
lowest  of  the  main  thermocline,  the  sound  speed  reaches  its  minimum  value;  which  can 
also  be  the  axis  of  the  sound  channel.  Figure  26  depicts  the  sound  channel.  According  to 
Snell’s  law,  as  sound  passes  through  different  water  layers,  the  ray  will  bend  toward  the 
region  where  sound  speed  is  slower,  thus  creating  a  propagation  path.  The  sound  speed 
reaches  its  minimum  at  the  axis  of  the  thermocline.  Therefore,  sound  will  bend  toward 
the  axis  of  the  thermocline. 
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Figure  25.  Three  Main  Layers  in  the  Ocean.  Adapted  from  Kalmbach  (2011). 


During  the  winter  seasons,  there  are  very  strong  winds  in  the  Mediterranean 
region.  Due  to  strong  winds  and  cooling,  a  large  scale  mixing  in  winter  makes  the 
temperature  of  the  surface  layer  almost  vertically  constant.  Under  the  seasonal  layer,  the 
sound  speed  can  decrease  rapidly  with  depth  within  the  main  thermocline.  Beneath  the 
main  thermocline,  a  deep  isothermal  layer  exists.  For  the  isothermal  case,  sound  speed 
starts  to  increase  due  to  pressure  effects.  These  cases  may  cause  a  shadow  zone,  which 
Figure  27  depicts. 
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Figure  26.  Sound  Channel.  Adapted  from  the  Office  of 
Ocean  Exploration  and  Research  (2001). 


Figure  27.  Shadow  Zone.  Adapted  from  DOSITS  (2014). 
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D.  SOUND  SPEED  PROFILES 


The  total  and  monthly  mean  SSPs  for  the  SMG-WOD  and  GDEM  are  calculated 
at  five  points  of  interest  (Figure  28)  selected  in  the  Western  (“E”),  Central  (“B”  and  “D”) 
and  Eastern  (“A”  and  “C”)  Mediterranean  Sea  to  cover  the  whole  basin  and  analyze  the 
regional  variations.  Table  14  provides  the  coordinates,  bottom  sediment  types  and  depths 
for  these  points.  Seasonal  variability  in  SMG-WOD  and  GDEM  SSPs  is  analyzed  for 
each  point. 


Figure  28.  Locations  of  Analysis  Points 


Table  14.  Points  A  through  E 


Location 

Latitude  (N) 

Longitude  (E) 

Depth  (m) 

Bottom 

Sediment 

Region 

A 

33.7794 

34.939 

1000 

clay 

Eastern 

B 

34.2399 

13.6275 

200 

very  fine  sand 

Central 

C 

32.1361 

27.9229 

1000 

sandy  mud 

Eastern 

D 

36 

18 

1000 

very  fine  silt 

Eastern 

E 

40.3946 

1,1739 

70 

sandy  mud 

Western 

1.  Point  A  (Eastern  Mediterranean  Sea) 

Summer  and  fall  seasons  exhibit  a  sharp  negative  gradient  in  sound  speed  due  to 
strong  surface  heating,  extending  down  to  approximately  80  meters.  Then,  the  sound 
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speed  continues  to  exhibit  a  weaker  negative  gradient  down  to  the  axis  of  the  sound 
channel,  layered  at  a  depth  of  400  meters  where  the  minimum  sound  speed  (1515  m/s) 
exists.  Beneath  the  sound  channel,  between  400  and  1000  m,  the  sound  speed  shows  a 
positive  gradient  due  to  increased  pressure. 

The  total  mean  for  SMG-WOD  SSPs  and  GDEM  SSPs  at  point  A  located  in  the 
Eastern  Mediterranean  Sea  are  plotted  in  Figure  29a.  The  corresponding  SMG-WOD 
multi-year  monthly  averages  and  standard  deviations  are  shown  in  Figure  30  along  with 
the  GDEM  monthly  mean  SSP. 
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Figure  29.  Total  Mean  SSP 


SMG-WOD  SSP  multi-year  monthly  average  and  GDEM  monthly  means 
correspond  well  in  surface  layers  and  isothermal  deep  layers.  Below  the  surface  layers  to 
the  depth  of  250  meters,  the  SMG-WOD  SS  is  about  3  m/s  higher  than  the  GDEM 
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monthly  mean.  The  SMG-WOD  standard  deviation  is  rather  small  as  expected.  However, 
it  does  not  mean  that  the  inter-annual  variability  is  small  since  the  standard  deviation  is 
not  a  full  representation  of  possible  outliers  or  non-periodic  inter-annual  and  decadal 
oscillations. 

January  through  March  exhibits  well-mixed  surface  layers  due  to  elevated  mixing. 
Surface  waters  are  not  heated  as  much  as  in  summer  and  fall.  Thus,  the  mean  sound  speed 
at  the  surface  is  lower  in  winter.  Below  the  surface  waters  to  the  depth  of  400  m,  the 
sound  speed  exhibits  a  negative  gradient.  Below  400  m,  the  sound  speed  shows  a  positive 
gradient  due  to  pressure  effects.  The  sound  channel  is  located  at  400  m  in  August. 
Convergence  zones  are  resent  during  all  seasons. 

The  greatest  variance  between  monthly  mean  and  total  mean  exists  in  the  upper 
layers,  particularly  from  January  through  March.  The  greatest  variance  is  about  12  m/s  in 
March  and  February.  It  is  about  9  m/s  in  the  January.  In  the  summer  and  fall  ,  there  is  a 
slight  difference  between  the  monthly  mean  and  total  mean  in  all  layers. 

2.  Point  B  (Central  Mediterranean  Sea) 

The  total  means  for  SMG-WOD  SSPs  and  GDEM  SSPs  at  point  B  located  in  the 
Central  Mediterranean  Sea  are  plotted  in  Figure  29b.  The  corresponding  SMG-WOD 
multi-year  monthly  averages  and  standard  deviations  are  shown  in  Figure  3 1  along  with 
the  GDEM  monthly  mean  SSP. 

January  through  April  exhibit  an  almost  uniform  gradient  from  the  surface  to  the 
bottom  layer.  This  was  assessed  as  slightly  positive  gradient.  Below  100  meters,  it  shows 
a  very  week  negative  gradient  due  to  temperature  effects.  In  April,  it  exhibits  a  very  weak 
negative  gradient  down  to  60  meters  resulting  from  some  surface  heating.  Beneath  60 
meters,  a  very  week  positive  gradient  exists.  The  surface  duct  is  present  from  December 
to  April  due  to  lower  sound  speed  at  the  surface.  In  addition,  we  expect  to  see  the  sound 
bending  downward  from  May  to  November  due  to  high  sound  speed  at  the  surface. 
Because  the  sea  is  shallow  there,  this  case  will  create  bottom  bounce  profiles. 
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Figure  30.  Mean  Sound-Speed  Profiles  (Point  A) 
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Figure  31.  Mean  Sound-Speed  Profiles  (Point  B) 


In  May,  surface  heating  starts  and  the  heat  increases  into  summer  and  early  fall. 
This  creates  a  very  strong  negative  gradient  in  the  upper  layer,  especially  in  August. 
Sound  speed  decreases  with  depth  rapidly  down  to  60  m.  Then  it  shows  a  very  weak 
negative  gradient  or  uniform  distribution  down  to  the  bottom. 

The  maximum  sound  speed  (1545m/s)  is  found  at  the  surface  in  the  August.  The 
minimum  sound  speed  (1508  m/s)  is  found  at  the  surface  in  February.  The  greatest 
variance  between  the  monthly  means  and  total  means  exist  in  the  upper  layers, 
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particularly  in  August  and  September  with  14  m/s  and  in  February  and  March  with  13 
m/s. 

Weak  surface  duct  propagation  can  be  expected  in  January  through  March  due  to 
slightly  positive  gradient.  The  rays  will  bend  upward  and  reflect  through  the  bottom.  This 
reflection  can  create  trapping  layers  above  100  m.  The  sound  can  propagate  over  a  long 
distance  in  trapping  layers.  This  is  a  very  disadvantageous  environment  for  submarines 
that  sail  at  periscope  depth.  Figure  32  illustrates  the  surface  duct  propagation. 


SURFACE  DUCT  PROPAGATION  PATH 


Figure  32.  Surface  Duct  Propagation.  Adapted  from  the  Department  of 
Oceanography,  Naval  Postgraduate  School  (2015). 


3.  Point  C  (Eastern  Mediterranean  Sea) 

The  total  means  for  SMG-WOD  SSPs  and  GDEM  SSPs  at  point  C  located  in  the 
Eastern  Mediterranean  Sea  are  plotted  in  Figure  29c.  The  corresponding  SMG-WOD 
multi-year  monthly  averages  and  standard  deviations  are  shown  in  Figure  33  along  with 
the  GDEM  monthly  mean  SSP. 

The  summer  season  exhibits  a  sharp  negative  gradient  in  the  surface  layer.  The 
sound  speed  decreases  rapidly  from  the  surface  to  80  m.  Then  a  thermocline  occurs 
below  80  m  extending  to  a  depth  of  400  m,  where  the  minimum  sound  speed  exists.  At 
the  lowest  depth  of  the  thermocline,  a  sound  channel  is  observed.  In  June,  the  sound 
speed  reaches  a  minimum  at  two  levels,  80  and  400  meters.  Thus,  we  expect  to  see  sound 
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channels  at  these  depths.  Having  two  sounds  channels  is  a  strong  indicator  that  sound  will 
propagate  over  a  long  distance. 

A  weak  positive  gradient  is  observed  in  the  upper  layer  in  January,  February  and 
March.  Positive  gradient  increases  from  January  to  March.  In  March,  the  positive 
gradient  extends  down  to  200  m. 

Below  the  upper  layer  with  the  positive  gradient  of  sound  speed,  a  negative 
gradient  is  observed.  The  minimum  sound  speed  is  reached  at  the  depth  of  400  meters.  In 
February  and  March,  there  are  two  minimum  sound  speed  levels,  at  the  surface  and  400 
meters.  Below  400  meters,  sound  speed  increases  constantly  with  depth.  This 
environment  may  cause  a  surface  duct  propagation  at  the  surface. 

The  greatest  variance  between  monthly  mean  total  means  exist  in  the  upper 
layers,  particularly  in  February  and  March.  Again,  there  is  a  difference  between  monthly 
means  in  the  thermocline,  which  must  be  the  result  of  different  vertical  resolutions 
between  the  two  datasets. 
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Figure  33.  Climatological  Monthly  Mean  Sound  Speed  Profiles  (Point  C) 
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4. 


Point  D  (Eastern  Mediterranean  Sea) 


The  total  means  for  SMG-WOD  SSPs  and  GDEM  SSPs  at  point  D  located  in  the 
Eastern  Mediterranean  Sea  are  plotted  in  Figure  29d.  The  corresponding  SMG-WOD 
multi-year  monthly  averages  and  standard  deviations  are  shown  in  Figure  34  along  with 
the  GDEM  monthly  mean  SSP. 

A  negative  gradient  in  the  surface  layer  exists  April  through  November,  followed 
by  a  uniform  gradient  in  the  surface  layer  extending  to  20  meters  in  December.  Then  it 
shows  a  negative  gradient  down  to  about  80  meters.  The  sound  speed  decreases  rapidly 
from  surface  to  the  depth  of  80  meters.  The  strongest  negative  gradient  occurs  in  August 
due  to  surface  heating.  The  sound  speed  reaches  its  minimum  at  a  depth  of  80  meters  in 
all  months  except  for  February.  This  pattern  may  create  a  surface  sound  channel.  The 
minimum  sound  speed  occurs  at  the  surface  in  February.  Point  D  shows  a  positive 
gradient  in  all  depths  as  exhibited  in  the  Arctic  Ocean.  Figure  35  presents  the  sound 
propagation  in  the  Arctic  Ocean. 
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Figure  34.  Sound  Propagation  in  the  Arctic  Ocean 
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In  March,  the  SMG-WOD  shows  a  minimum  sound  speed  at  the  surface,  but  the 
GDEM  shows  it  at  a  depth  of  80  meters  as  in  other  months. 

The  greatest  variance  (approximately  17  m/s)  between  monthly  and  total  means 
exists  in  the  upper  layers,  particularly  from  January  to  March. 
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Figure  35.  Climatological  Monthly  Mean  Sound  Speed  Profiles  (Point  D) 
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5.  Point  E  (Western  Mediterranean  Sea) 

The  total  means  for  SMG-WOD  SSPs  and  GDEM  SSPs  at  point  E  located  in  the 
Western  Mediterranean  Sea  are  plotted  in  Figure  29e.  The  corresponding  SMG-WOD 
multi-year  monthly  averages  and  standard  deviations  are  shown  in  Figure  36  along  with 
the  GDEM  monthly  mean  SSP. 
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Figure  36.  Climatological  Monthly  Mean  Sound  Speed  Profiles  (Point  E) 
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In  December  through  March,  monthly  means  are  almost  uniform  due  to  strong 
and  cold  winds.  Wind  stirring  mixes  the  surface  layer  and  creates  an  isothermal  surface 
layer  of  water.  It  is  well-mixed  and  very  shallow.  Thus,  thermocline  does  not  exist  in 
these  months.  Surface  heating  starts  in  April  and  increases  until  August  as  expected. 
Maximum  mean  values  (1535  m/s)  are  found  in  August.  The  August  profile  shows  a 
stronger  negative  gradient.  The  SMG-WOD  standard  deviation  is  very  small  here,  and  the 
annual  SSPs  is  close  to  mean  monthly  values.  The  greatest  variance  (approximately  15 
m/s)  between  monthly  mean  and  total  mean  occurs  in  August  and  September  in  the  upper 
of  a  depth  of  30  meters.  The  minimum  variance  occurs  in  May  in  the  upper  10  m  layer 
and  in  June  the  deeper  layer  at  the  depth  of  20  m.  The  maximum  variance  is  observed  two 
months  later,  in  August,  due  to  significant  heating. 

Table  15  summarizes  expected  acoustic  propagation  types  and  months  when 
maximum  seasonal  variations  are  observed. 


Table  15.  Possible  Propagation  Types  and  Maximum  Variations 


Location 

Region 

Depth 

(m) 

Month 

Possible  Propagation  Type 

Maximum 

Variation 

A 

Eastern 

1000 

August 

Sound  Channel  (400  m) 

February  and 
March 

January 

Sound  Channel  (400  m) 

B 

Central 

200 

August 

Sound  Bends  Downward 

August 

February 

Weak  Surface  Duct 

C 

Eastern 

1000 

August 

Sound  Channel  (  400  m) 

February  and 
March 

March 

Surface  Duct 

D 

Eastern 

1000 

August 

Sound  Channel  (80  m) 

January  - 
March 

February 

Surface  duct 

E 

Western 

70 

August 

No  Sound  Channel  or  Duct 

August 

January 

No  Sound  Channel  or  Duct 
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At  points  A,  C  and  D  it  is  possible  to  see  convergence  zones  in  layers  deeper  than 
those  shown  in  Table  15.  Maximum  variance  in  the  Eastern  Mediterranean  Sea  occurs  in 
January  through  March;  for  the  central  and  the  Western  Mediterranean  Sea,  the 
maximum  variance  occurs  in  the  August. 

E.  MEAN  TRANSMISSION  LOSS  -GDEM 

Transmission  loss  has  been  modeled  using  Bellhop  for  a  3500  Hz  sound  source 
located  at  40  m  and  101  evenly  distributed  receiver  depths,  along  a  west-east  path  for 
maximum  propagation  range  of  70  nm,  using  GDEM  monthly  mean  sound  speed  profiles. 
To  be  able  to  see  how  sound  propagates  in  the  Mediterranean  Sea  will  provide  us  with 
information  on  where  the  submarines  can  be  detected.  Figures  37-41  show  transmission 
loss  for  each  month  at  points  A-E. 

a.  Point  A 

We  expected  to  see  a  sound  channel  at  a  depth  of  400  m  in  January  and  August; 
however,  Figure  37  shows  no  sound  channel  in  August.  In  August,  sound  rays  propagated 
downward  to  1600  m  due  to  a  very  high  sound  speed  at  the  surface.  Then,  they  refracted 
upward  due  to  high  speed  at  a  depth  of  1600  meters.  Therefore,  sound  rays  created 
convergence  zone  at  the  surface  at  an  approximately  35  km  range  and  submarines  could 
have  been  detected  from  long  distances  where  the  convergence  path  occurred.  Also,  there 
was  a  direct  path  near  the  surface  within  an  8  km  range,  and  the  TL  was  about  80-85  dB. 
Thus,  submarines  could  also  have  been  detected  near  the  surface  within  an  8  km  range. 
Between  8  km  and  35km,  it  is  difficult  to  detect  submarines  due  to  the  very  large  shadow 
zone. 

For  January,  we  detected  a  surface  duct  propagation  within  the  17  km  range  and  a 
convergence  zone  at  the  35  km  range.  Thus,  submarines  could  have  been  detected  within 
the  17  km  range  at  the  surface.  They  also  could  have  been  detected  from  long  distances 
where  the  convergence  path  occurred.  For  February  and  March,  we  detected  a  surface 
duct  within  the  45  km  range  and  a  sound  channel  convergence  zone  at  a  depth  of  750- 
1700  meters.  The  convergence  path  occurred  at  the  32  km  range.  Thus,  submarines  could 

have  been  detected  at  the  surface  from  long  distances. 
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Figure  37.  Monthly  Mean  TL  Calculated  from  BELLHOP  model  at  Point  A 
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Figure  37  (cont’d.).  Monthly  Mean  TL  Calculated  from  BELLHOP  model  at  Point  A 
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For  September  and  October,  bottom  bounce  occurred.  Sound  struck  the  bottom 
and  reflected  upward  to  the  surface  at  the  25-32  km  range.  For  June  and  July,  there  was  a 
very  large  shadow  zone,  which  extended  to  a  depth  of  700  meters  within  the  32  km  range. 
These  results  mostly  correspond  to  the  convergence  zone  expected  from  mean  sound 
speed  profiles,  but  not  for  sound  channel 

These  convergence  zones  are  of  main  concern  for  USW.  There  are  many 
submarine  detection  tactics  based  on  convergence  zones.  Point  A  provides  a  high 
possibility  of  detecting  submarines  at  long  distances. 

b.  Point  B 

We  expected  to  see  a  surface  duct  in  February  and  sound  to  bend  downward  in 
August.  From  January  through  April,  a  surface  duct  propagation  was  modeled  along  the 
long  distances  (Figure  38).  After  April,  the  sound  started  to  bend  downward  due  to 
heating  at  the  surface  and  bottom-bounce  occurred  due  to  the  shallow  water  environment. 
December  exhibited  both  a  surface  duct  and  bottom  bounce.  These  results  agrees  with 
what  expected  from  mean  sound  speed  profiles. 
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Figure  38.  Same  as  Figure  37  except  for  Point  B 
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Figure  38  (cont’d.). 


Same  as  Figure  37  except  for  Point  B 
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Figure  38  (cont’d.).  Same  as  Figure  37  except  for  Point  B 


c.  Point  C 

We  expected  to  see  surface  duct  in  March  and  a  sound  channel  at  400  meters  in 
August.  We  also  expected  the  convergence  zone  due  to  deep  layer  in  Point  C.  From 
December  through  April,  there  are  a  surface  duct  and  convergence  zones  in  each  figure 
(Figure  39).  The  strongest  surface  duct  occurred  in  March  as  we  expected.  However,  we 
did  not  see  a  sound  channel  in  August  as  we  expected.  From  May  through  November,  a 
convergence  zone  but  no  sound  channel  was  detected.  The  upper  limiting  ray  of 
convergence  zones  varied  between  500  and  1700  meters;  maximum  upper  limiting  ray  of 
the  convergence  zone  took  place  in  September  and  October,  and  the  minimum  upper 


66 


limiting  ray  of  the  convergence  zone  occurred  in  March.  Shadow  zones,  created  by 
convergence  zones,  would  have  been  a  good  place  to  avoid  detection. 
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Figure  39.  Same  as  Figure  37  except  for  Point  C 
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Figure  39  (cont’d).  Same  as  Figure  37  except  for  Point  C 
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Figure  40.  Same  as  Figure  37  except  for  Point  D 
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Figure  40  (cont’d.).  Same  as  Figure  37  except  for  Point  D 
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From  January  through  March,  a  surface  duct  and  a  deep  convergence  zone 
occurred.  These  results  agree  with  our  expectation  from  the  mean  sound  speed  profiles  as 
the  Arctic  Ocean  profile  exhibits  surface  duct  and  convergence  zones  together  as  long  as 
there  are  deep  layers.  From  April  through  December,  we  expected  to  see  a  sound  channel 
at  a  depth  of  80  meters  in  conjunction  with  convergence  zones.  However,  because  of  a 
very  sharp  negative  gradient  from  surface  down  to  80  meters,  where  sound  speed  reaches 
its  minimum,  we  only  observed  convergence  zones.  The  upper  limiting  ray  of 
convergence  zones  varied  between  300  and  1200  meters;  the  maximum  upper  limiting 
ray  of  convergence  zone  took  place  in  September  and  October,  and  the  minimum  upper 
limiting  ray  of  the  convergence  zone  occurred  in  April. 


e.  Point  E 

In  the  shallow  water,  a  direct  path  was  modeled  (Figure  41).  Sound  bends 
downward  and  bounced  from  the  bottom  to  the  surface.  The  20  km  range  exhibited  a 
deep  layer,  and  a  surface  duct  emerged  from  January  through  March  due  to  low  speed  at 
the  surface.  In  April,  surface  sound  channels  took  place  at  the  100  meters.  From  May 
through  December,  the  sound  channel  deepens,  and  a  weak  sound  channel  occurs  at  200 
meters.  However,  the  TL  was  significant  in  these  sound  channels.  The  surface  duct  and 
sound  channels  are  not  as  strong  as  those  occurring  in  January  through  May.  Thus,  would 
have  been  possible  to  detect  submarines  at  the  surface  at  very  long  distances  from 
January  through  April,  but  not  from  May  through  December. 
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Figure  41 .  Same  as  Figure  37  except  for  Point  E 
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Figure  41  (cont’d.). 


Same  as  Figure  37  except  for  Point  E 
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Figure  41  (cont’d.).  Same  as  Figure  37  except  for  Point  E 
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2. 


GDEM  AND  SMG-WOD  COMPARISON 


January  and  August  have  been  chosen  to  analyze  the  range  and  patterns  of  sound 
propagation,  as  well  as  TL.  In  Figures  41-46,  GDEM  modeling  results  are  shown  in  the 
left  panels,  and  the  SMG-WOD  modeling  results  are  shown  in  the  right  panels. 
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Figure  42.  Point  A  January  (upper)  and  August  (lower) 


At  the  point  A  location  in  January  (Figure  42,  upper  panels),  the  sound 

propagation  patterns  were  similar  to  each  other.  The  SMG-WOD  (right)  exhibited  surface 

duct  propagation  along  the  45  km  range.  However,  the  GDEM  (left)  only  showed  a 

surface  duct  up  to  the  17  km  range,  but  it  exhibited  another  surface  duct  after  the 

convergence  path  at  the  32  km  range.  The  upper  and  lower  limiting  rays  of  convergence 
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between  the  datasets  were  in  agreement.  However,  there  was  a  slight  difference  between 
their  TLs  in  the  shadow  zones.  The  GDEM  reflects  a  slightly  lower  TL  in  the  shadow 
zones.  In  August  (Figure  42,  lower  panels),  sound  propagation  patterns  were  similar  to 
each  other.  However,  there  was  a  slight  difference  in  TL  in  the  direct  path.  The  upper 
limiting  ray  of  the  convergence  zone  was  150  meters  deeper  in  SMG-WOD  (right)  than 
in  the  GDEM  (left).  Lower  TL  is  modeled  with  SMG-WOD  in  shadow  zones.  In  August 
(Figure  43,  lower  panels)  higher  attenuation  is  modeled  at  surface,  and  lower  attenuation 
is  modeled  in  the  deep  layer.  Both  the  datasets  produced  shorter  propagation  ranges  in 
August  than  in  January. 
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Figure  43.  Point  B  January  (upper)  and  August  (lower) 


75 


At  the  point  B  location  in  January  (Figure  43,  upper  panels),  sound  propagation 
patterns  were  similar  to  each  other.  Both  models  produced  a  surface  duct  at  the  same 
depth  and  range.  After  the  22  km,  a  bottom  bounce  pattern  was  observed  for  both  the 
models;  however,  GDEM  (left)  exhibited  a  clearer  bottom  bounce  regime.  TL  was  almost 
the  same  in  both  the  models.  Additionally,  SMG-WOD  (right)  exhibited  a  surface  duct 
with  thicker  profiles.  In  August  (Figure  43,  lower  panels),  the  sound  propagation  patterns 
were  very  much  like  for  both  the  datasets.  The  bottom  bounce  pattern  was  observed  in 
each  figure  and  there  was  no  surface  duct  after  the  38  km  range  at  the  surface.  The 
GDEM  (left)  exhibited  a  clearer  bottom  bounce  regime,  and  TL  was  somewhat  higher  in 
GDEM.  The  SMG-WOD  bottom  bounce  pattern  (right)  propagates  further  than  the 
GDEM  with  thicker  sound  profiles.  In  August,  higher  attenuation  values  occurred  at  the 
surface  and  lower  attenuation  occurred  in  the  deep  layer,  which  produces  shorter  ranges 
than  in  January. 


500 

1000 

E 

^"1500 

O- 

© 

°  2000 
2500 

0000 


BELLHOP-  Ray090 
Freq  =  3500  Hz  Sd  =  40  m 


2  4  6  8  10 

Range  (m) 

BELLHOP-  Ray090 
Freq  =  3500  Hz  Sd  =  40  m 


2  4  6  8  10 

Range  (m) 


12 

*104 


12 

xlOJ 


50 

55 

50 

65 

70 

75 

80 

85 

90 

95 

100 


500 


1000 


1500 


2000 


2500 


3000 


12 

xIO4 


12 

xIO4 


500 


1000 


31500 


®  2000 


2500 


3000 


2  4  6  8  10 

Range  (m) 


500 


1000 


3 1500 

s 

§2000 


2500 


3000 


BELLHOP-  Ray090 
Freq  =  3500  Hz  Sd  =  40  m 


2  4  6  8  10 

Range  (m) 

BELLHOP-  Ray090 
Freq  =  3500  Hz  Sd  =  40  m 


Figure  44.  Point  C  (upper)  and  August  (lower) 
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At  the  point  C  location  in  January  (Figure  44,  upper  panels),  both  the  datasets 
produced  the  similar  surface  duct  and  convergence  zones  with  only  a  slight  difference  in 
TL  (exaggerated  because  of  different  colorbars).  Additionally,  the  SMG-WOD  (right) 
showed  a  surface  duct  with  thicker  profiles.  In  August,  both  the  models  produced  similar 
convergence  zones  with  a  light  difference  in  TL.  In  August,  there  are  higher  attenuation 
values  at  the  surface  but  lower  attenuation  in  the  deep  layer,  which  results  in  shorter 
ranges  than  January. 
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Figure  45.  Point  D  (upper)  and  August  (lower) 


At  the  point  D  location  in  January  (Figure  45,  upper  panels),  both  the  datasets 
produced  a  surface  duct,  sound  channel  and  convergence  zones;  however,  there  are  some 
notable  differences.  The  maximum  lower  limiting  ray  of  the  convergence  zone  extended 
deeper  in  the  SMG-WOD  (right)  at  the  22  km  range.  There  was  also  a  higher  the  TL  in 

shadow  zones  in  GDEM  (left).  After  the  40  km  range,  there  were  different  types  of 
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convergence  zones  in  the  SMG-WOD.  In  August  (Figure  45,  lower  panels),  the 
maximum  lower  limiting  ray  of  the  convergence  zone  extended  deeper  in  the  SMG-WOD 
at  the  22  km  range.  The  GDEM  had  a  higher  TL  in  deep  shadow  zones  but  lower  TL  in 
surface  shadow  zones.  The  maximum  lower  limiting  ray  of  the  convergence  zone 
extended  deeper  in  SMG-WOD  than  in  the  GDEM  in  the  bottom  layer 
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Figure  46.  Point  E  (upper)  and  August  (lower) 


At  the  point  E  location  in  August  (Figure  46,  lower  panels),  the  sound 
propagation  patterns  were  similar  to  each  other.  Both  of  the  figures  show  a  surface  sound 
channel  and  convergence  zone.  The  direct  path  took  place  in  shallow  water.  Then,  a  weak 
sound  channel  occurred  at  200  meters.  The  SMG-WOD  (right)  exhibited  lower  TL  in  the 
sound  channel  and  higher  TL  in  the  convergence  zone.  The  lower  limiting  ray  of 
convergence  zones  penetrated  deeper  in  GDEM  than  in  SMG-WOD. 
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VI.  INTERANNUAL  VARIABILITY 


A.  TIME  VARIATIONS  OF  (T,  S,  SS)  PROFILES 

Temperature,  salinity,  and  sound  speed  profiles  at  the  five  selected  locations  are 
described  in  this  section.  The  (T,  S)  data  are  obtained  from  the  SMG-WOD  dataset  and 
the  SS  profiles  are  calculated  using  the  Medwin  equation. 

1.  Temperature 

Figures  47-51  showed  temporal  variability  of  SMG-WOD  temperature  from 
1960-2014  with  one-month  resolution  for  points  A-E.  The  temperature  profiles  at  all  the 
locations  exhibit  strong  seasonality  associated  with  surface  heat  fluxes  and  variations  in 
the  vertical  structure  of  water  column.  However,  this  seasonality  was  strongly  modulated 
by  inter-annual  and  decadal  variability  as  described  below  for  each  location. 


I960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010  2015 

Figure  47.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  A 
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At  point  A  (Figure  47)  the  inter-annual  variability  is  evident  in  the  variations  of 
summer  temperatures  at  sea  surface,  the  depth  of  the  mixed  layer  and  the  strength  of  the 
thermocline.  The  temperature  was  colder  between  1980  and  1990  at  the  surface  and  in  the 
intermediate  layer,  with  shallower  summer  mixed  layer  observed  in  1985-1986.  The 
temperature  increased  between  1990  and  2014.  Temperature  between  1960  and  1980  was 
also  warmer  than  1980  to  1990,  however,  it  was  not  warmer  than  period  between  1990 
and  2014.  Strong  variations  of  temperature  with  higher  temperatures  propagating  down  to 
500  meters  were  observed  in  2009-2010. 


Temperature  at  Point  B 


200 


1960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010  2015 


30 

28 

26 

24 

22 

20 


18 


16 

14 


Figure  48.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  B 

Point  B  shows  some  important  variations  between  specific  years  (Figure  48). 
Temperature  at  the  surface  between  1997  and  1998  was  warmer  than  in  other  years.  Also, 
warmer  water  between  2007  and  2008  extended  down  to  180  meters. 
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Figure  49.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  C 

At  Point  C  (Figure  49)  the  temperature  at  the  surface  exhibited  a  strong  decadal 
variability.  Especially  between  2000  and  2014,  the  surface  water  gets  warmer  and  this 
warmer  water  penetrate  to  deeper  layers  than  in  other  years.  An  anomalously  deep 
penetration  of  warm  water  down  to  300  meters  was  observed  during  summer  1992. 
Interestingly,  it  was  accompanied  by  colder  than  usual  surface  temperature  and  rather 
shallow  mixed  layer  at  the  end  of  spring.  An  anomalously  warm  surface  water  but  no 
deepening  of  the  warm  surface  layer  was  observed  during  the  summer  2012 
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Figure  50.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  D 

At  Point  D  (Figure  50),  a  strong  contrast  existed  in  summer  surface  temperatures 
before  and  after  1980.  The  temperature  at  the  surface  started  to  increase  after  1980. 
Before  the  1980,  the  surface  water  is  about  24  °C.  Then  it  gets  warmer  at  the  surface  and 
reaches  29  °C.  Between  2000  and  2014  years,  the  warmer  water  starts  to  penetrate  down 
to  400  meters. 
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Figure  5 1 .  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  E 

Temperature  increased  between  1990  and  2014.  The  waters  between  surface  and  a 
depth  of  20  meters  reached  the  higher  temperature  than  previous  years;  however,  the 
warmer  water  penetrated  to  deeper  layer  between  1997  and  1998.  Sea  surface 
temperatures  between  1960  and  1990  were  lower  than  later  decades. 

2.  Salinity 

Figures  52-56  show  temporal  variability  of  SMG-WOD  salinity  from  1960-2014 
with  one-month  resolution  for  points  A-E.  While  some  seasonality  is  seen  in  the  salinity 
profiles,  the  inter-annual  and  decadal  variability  is  much  stronger,  most  likely  associated 
with  inter-annual  and  decadal  variability  of  freshwater  fluxes  and  water  mass  formation. 


83 


Salinity  at  Point  A 

50 
100 

150 

200 
E 

^250 
0) 

Q 

300 
350 
400 
450 
500 

1960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010  2015 

Figure  52.  Temporal  Variability  of  SMG-WOD  Salinity  at  Point  A 

The  salinity  showed  inter-annual  variations  between  1960  and  2014  (Figure  52). 
Overall,  point  A  is  becoming  saltier  from  1960s.  The  inter-annual  and  decadal  variability 
at  this  location  seems  to  be  synchronized  with  the  inter-annual  variability  of  temperature 
(Figure  47).  Specifically,  the  salinity  increased  between  1990  and  1995,  and  between 
2008  and  2014.  It  did  not  increase  between  1995  and  2008. 
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Figure  53.  Temporal  Variability  of  SMG-WOD  Salinity  at  Point  B 
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Point  B  has  saltier  water  below  about  80  m  (Figure  53).  That  must  be  due  to  less 
saline  Modified  Atlantic  Water  at  the  surface.  The  salty  layer  exhibited  a  very  robust 
seasonality  with  an  episodic  increase  in  salinity  in  2007.  At  the  surface,  there  is  high 
inter-annual  variations  in  salinity;  however,  the  salinity  did  not  show  steady  increase 
between  1960  and  2014.  The  salinity  decreased  at  the  surface  at  some  years  between 
1970  and  1990.  These  values  were  lower  than  values  between  1960  and  1970.  After 
1990,  the  salinity  increased,  especially  between  2007  and  2008. 
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Figure  54.  Inter-annual  Salinity  Variations  at  Point  C 

At  point  C,  a  three-layer  structure  was  observed  with  highly  seasonal  upper  layer 
down  to  about  50  meters  on  average,  salinity  minimum  between  50  and  100  meters  and  a 
saltier  water  below  (Figure  54).  The  upper  layer  demonstrated  higher  salinity  during 
summer  associated  with  high  evaporation.  The  inter-annual  variability  in  surface  heat  and 
freshwater  fluxes  drove  inter-annual  variability  in  the  upper  layer,  except  between  1990 
and  2005.  2005-2010  was  probably  the  most  anomalous  period  with  extremely  high 
salinity  at  sea  surface  and  several  fresh  water  intrusions  below  50  meters.  After  2005, 
salinity  also  started  to  increase  below  150  meters. 
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Figure  55.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  D 

At  point  D  (Figure  55),  there  is  an  apparent  two-layer  haline  structure  with  fresher 
layer  of  about  100  meters  thick  laying  on  top  of  the  water  column.  The  two  layers  were 
separated  by  rather  strong  halocline.  While  some  seasonality  is  evident  in  the  upper  layer, 
there  were  high  inter-annual  variations  between  surface  100-meter  layer  and  the  rest  of 
the  water  column.  The  salinity  at  the  surface  decreased  between  1980  and  1995.  After 
1995,  the  salinity  started  to  increase.  It  reached  its  maximum  value  at  2007.  The  salinity 
increases  in  all  layers  between  2010  and  2014.  The  deep  layer  showed  its  highest  inter¬ 
annual  variations  between  2010  and  2014.  In  2005  an  intrusion  of  salty  water  was  seen  in 
the  subsurface  layer  between  30  and  200  meters. 
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Figure  56.  Temporal  Variability  of  SMG-WOD  Temperature  at  Point  E 

Point  E  (Figure  56)  exhibits  high  inter-annual  variability  in  all  depths.  After  1975, 
the  salinity  started  to  increase  until  1990.  It  decreases  rapidly  between  1990  and  1995. 
After  the  1995,  the  salinity  started  to  increase  again  in  all  layers.  The  salinity  reached  its 
maximum  value  between  2010  and  2014.  Besides  these  long-term  trends,  occasional 
intrusions  of  anomalously  fresh  water  in  1980,  1993  and  2011. 

3.  Sound  Speed 

Figures  57-61  show  temporal  variability  of  SMG-WOD  sound  speed  from  1960- 
2014  with  one-month  resolution  for  points  A-E.  Comparison  of  sound  speed  (Figures  57- 
61)  and  temperature  (47-51)  variability  shows  that  the  upper  layer  sound  speed  is 
determined  by  temperature  as  expected,  while  the  lower  level  is  determined  by  the 
pressure  and  shows  neither  seasonal  nor  inter-annual  variability.  It  is  however  expected 
that  strong  inter-annual  variability  of  salinity  will  affect  the  inter-annual  variability  of 
sound  speed. 
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Figure  57.  Temporal  Variability  of  SMG-WOD  Sound  Speed  at  Point  A 


At  point  A  (Figure  57)  the  sound  speed  at  the  surface  was  almost  same.  Below  a 
depth  of  40  meters,  the  sound  speed  shows  higher  inter-annual  variations.  The  sound 
speed  decreases  between  1980  and  1990.  The  sound  speed  reaches  its  maximum  values 
between  2000  and  2014.  These  results  agreed  with  temperature  and  salinity  results  of 
Point  A  because  both  of  them  increases  between  2000  and  2014. 
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Figure  58.  Temporal  Variability  of  SMG-WOD  Sound  Speed  at  Point  B 

The  Point  B  (Figure  58)  show  some  higher  variations  in  temperature  and  salinity 
in  specific  years.  Higher  sound  speed  penetrated  to  a  deeper  layer  between  2007  and 
2008.  Sound  speed  reached  its  maximum  at  the  surface  between  1995  and  1998. 
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Figure  59.  Temporal  Variability  of  SMG-WOD  Sound  Speed  at  Point  C 

At  point  C  (Figure  59)  sound  speed  shows  decadal  oscillations  at  the  surface. 
There  were  higher  variations  below  80  meters.  After  2005,  we  found  higher  sound  speed 
at  this  depth.  Also,  we  found  higher  sound  speed  at  the  same  depth  in  1992.  The  sound 
speed  also  increased  at  the  surface  between  2000  and  2014.  The  sound  speed  reached  its 
minimum  value  between  1980  and  1985. 
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Figure  60.  Temporal  Variability  of  SMG-WOD  Sound  Speed  at  Point  D 


At  point  D  (Figure  60)  the  sound  speed  decreased  between  1980  and  2000,  with 
an  exception  of  1990  and  1992.  There  are  very  low  sound  speed  values  at  deeper  layers 
during  these  years.  From  2000  to  2014,  the  sound  speed  increased  at  all  depths,  especially 
between  60  and  200  meters. 
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Figure  6 1 .  Temporal  Variability  of  SMG-WOD  Sound  Speed  at  Point  E 


At  point  E  (Figure  61)  the  sound  speed  between  1960  and  1990  is  lower  than 
during  later  decades.  It  starts  to  increase  slightly  after  1990  and  reaches  its  maximum 
between  2010  and  2015. 

B.  EOF  ANALYSIS 

Empirical  orthogonal  functions  (EOF)  analysis  was  conducted  to  identify  the 
inter-annual  variability  of  the  sound  speed  profiles  relative  to  the  multiyear  averaged 
profiles,  c(z,.,f  ),  from  the  SMG-WOD  data  between  1960  to  2014  are  expanded  into 

EOFs, 

c{zt  JP)  =  X  PCi  (?p)Ej  (*« )  > 


93 


where  £j(z,)  is  the  jth  EOF  mode  with  the  unit  of  m/s;  and  PC  .(?  )  is  the  jth  principal 

component  without  unit.  Here,  j  =  1,  2,  J;  J  is  truncated  number  of  the  EOF  modes. 
Table  16  shows  variance  contained  in  the  four  modes  at  selected  locations. 


Table  16.  Variances  of  First  Four  EOF  Modes 


Point 

Mode  1 

%  Variance 

Mode  2 
%  Variance 

Mode  3 
%  Variance 

Mode  4 
%  Variance 

Cumulative 

Variance 

A 

34.1 

18.9 

13.8 

6.1 

72.9 

B 

39.9 

17.7 

9.4 

8.3 

75.3 

C 

42.2 

18.3 

9.6 

7.1 

77.2 

D 

43.9 

26.3 

9.8 

4.1 

84.1 

E 

67.2 

19.8 

5.3 

3.8 

96.1 

Figures  62  to  66  show  the  eigenvalues  and  cumulative  variance  as  percentage  of 
total  variability  from  the  first  to  truncated  (30st)  mode  for  the  five  locations  (A,  B,  C,  D, 
E).  In  particular,  at  the  point  A,  it  needs  10-1 1  modes  to  reach  90%  of  total  variance.  At  a 
shallower  region  (point-E),  three  modes  are  required  to  reach  90%  of  variance.  Point  D, 
which  is  located  to  the  west  of  Point  A  and  at  the  same  bottom  depth  as  Point  C,  demands 
six  modes  to  reach  90%  of  the  total  variance.  However,  this  result  may  change  according 
to  the  data  scarcity  and  depths  in  other  regions. 
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Figure  62.  Point  A  (Eastern):  Eigenvalues  and  Mode  Variance 
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Figure  63.  Point  B  (Central):  Eigenvalues  and  Mode  Variance 


95 


200 


Sound  Speed  &  A  =  o  =  321 361 


10  15  20 

Mode  Number 


Figure  64.  Point  C  (Eastern):  Eigenvalues  and  Mode  Variance 
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Figure  65.  Point  D  (Eastern):  Eigenvalues  and  Mode  Variance 
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Figure  67  and  68  shows  four  first  EOFs  and  corresponding  EOF  coefficients 
(Principal  Components)  for  Point  A.  Only  the  first  four  EOFs  are  shown  because  they 
explain  most  of  the  variability  in  the  data.  The  modes  contribute  with  different  sign. 
Table  17  through  21  demonstrate  the  results. 
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Figure  67.  The  First  Four  EOFs  of  Point  A 
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Figure  68.  EOF  Coefficients  of  Sound  Speed  of  Point  A 


According  to  the  EOF  analysis  of  Point  A,  we  found  that  sound  speed  does  not 
change  below  a  depth  of  800  meters  at  that  location.  We  also  found  that  there  was  a  very 
high  variability  from  1980  to  1990  and  from  2000  to  2014.  Therefore,  we  found  that 
sound  speed  at  the  surface  was  faster  between  1960  and  1980;  and  slower  between  1980 
and  1990. Those  two  periods  need  to  be  examined  in  much  greater  detail  to  determine 
what  dynamic  processes  were  behind  these  changes. 
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Table  17.  EOF  Analysis  at  Point  A 


Mode  Number 

EOF  Mode 

PC 

Sound  Speed 

Depth 

Value 

Time 

Value 

1 

20 

Positive 

1960-1980 

2000-2014 

Positive 

Increase 

80 

Positive 

Increase 

300 

Positive 

Increase 

800 

Zero 

No  change 

1 

20 

Positive 

1980-1990 

Negative 

Decrease 

80 

Positive 

Decrease 

300 

Positive 

Decrease 

800 

Zero 

No  change 

2 

20 

Positive 

1960-1980 

2000-2008 

Positive 

Increase 

80 

Positive 

Increase 

300 

Negative 

Decrease 

800 

Zero 

No  change 

2 

20 

Positive 

1980-1990 

2008-2010 

Negative 

Decrease 

80 

Positive 

Decrease 

300 

Negative 

Increase 

800 

Zero 

No  change 

3 

20 

Negative 

1960-1970 

Positive 

Decrease 

80 

Positive 

Increase 

300-500 

Negative 

Decrease 

800 

Zero 

No  change 

3 

20 

Negative 

2009-2010 

Negative 

Increase 

80 

Positive 

Decrease 

300-500 

Negative 

Increase 

800 

Zero 

No  change 

4 

20 

Positive 

1960-1970 

1980-1985 

Negative 

Decrease 

80 

Negative 

Increase 

300-500 

Negative 

Increase 

800 

Zero 

No  change 

4 

20 

Positive 

2009-2010 

Positive 

Increase 

80 

Negative 

Decrease 

300-500 

Negative 

Decrease 

800 

Zero 

No  change 
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Table  18.  EOF  Analysis  at  Point  B 


Mode  Number 

EOF  Mode 

PC 

Sound  Speed 

Depth 

Value 

Time 

Value 

1 

10 

Negative 

1960-1970 

1980-1990 

1997-1998 

Positive 

Decrease 

60 

Negative 

Decrease 

80 

Negative 

Decrease 

140 

Negative 

Decrease 

1 

10 

Positive 

2007-2008 

1970-1980 

Negative 

Decrease 

60 

Positive 

Decrease 

80 

Positive 

Decrease 

140 

Positive 

Decrease 

2 

10 

Negative 

1960-1980 

1997-1998 

2007-2008 

Positive 

Decrease 

60 

Positive 

Increase 

80 

Positive 

Increase 

140 

Positive 

Increase 

3 

10 

Negative 

1970-1980 

Positive 

Decrease 

60 

Positive 

Increase 

80 

Negative 

Decrease 

140 

Negative 

Decrease 

3 

10 

Negative 

1997-1998 

Negative 

Increase 

60 

Positive 

Decrease 

80 

Negative 

Increase 

140 

Negative 

Increase 

4 

10 

Positive 

1960-1970 

1980-1992 

1997-1998 

2007-2008 

Positive 

Increase 

60 

Positive 

Increase 

80 

Positive 

Increase 

140 

Negative 

Decrease 

According  to  the  EOF  analysis,  Point  B  has  a  very  high  variability  for  almost  all 
of  the  period  (Figures  69  and  70,  and  Table  18).  However,  2007-2008  and  1997-1998 
exhibited  the  greatest  variability.  Sound  speed  throughout  all  depths  decreased  during 
these  periods.  In  addition,  sound  speed  decreased  from  1960  to  1980.  Those  periods  need 
to  be  examined  further  to  determine  what  dynamic  processes  were  behind  these  changes. 
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Figure  7 1 .  The  First  Four  EOFs  of  Point  C 
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Figure  72.  Coefficients  of  Sound  Speed  of  Point  C 
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Table  19.  EOF  Analysis  at  Point  C 


Mode  Number 

EOF  Mode 

PC 

Sound  Speed 

Depth 

Value 

Time 

Value 

1 

20 

Positive 

1980-1995 

Negative 

Decrease 

80 

Positive 

Decrease 

300 

Positive 

Decrease 

600 

Zero 

No  change 

1 

20 

Positive 

2000-2014 

Positive 

Increase 

80 

Positive 

Increase 

300 

Positive 

Increase 

600 

Zero 

No  change 

2 

20 

Negative 

1968-1989 

2010-2014 

Positive 

Decrease 

80 

Negative 

Decrease 

300 

Positive 

Increase 

600 

Zero 

No  change 

2 

20 

Negative 

2000-2010 

Negative 

Increase 

80 

Negative 

Increase 

300 

Positive 

Decrease 

600 

Zero 

No  change 

3 

20 

Positive 

1965-1985 

2010-2014 

Positive 

Increase 

80 

Negative 

Decrease 

300 

Positive 

Increase 

600 

Zero 

No  change 

3 

20 

Positive 

2000-2010 

Negative 

Decrease 

80 

Negative 

Increase 

300 

Positive 

Decrease 

600 

Zero 

No  change 

4 

20 

Positive 

1982-2005 

Positive 

Increase 

80 

Negative 

Decrease 

300 

Negative 

Decrease 

600 

Zero 

No  change 

4 

20 

Positive 

2010-2014 

Negative 

Decrease 

80 

Negative 

Increase 

300 

Negative 

Increase 

600 

Zero 

No  change 

We  found  there  was  a  very  high  variability  from  1980  to  1995  and  from  2000  to 
2014  at  the  point  C  location  ((Figures  71  and  72,  and  Table  19).  From  1980  to  1995,  the 
sound  speed  in  the  upper  and  intermediate  layers  decreased  dramatically;  conversely,  the 
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upper  and  intermediate  layers  increased  from  2000  to  2014.  We  found  that  sound  speed 
does  not  change  below  a  depth  of  600  meters  during  the  entire  period. 
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Figure  73.  The  First  Four  EOFs  of  Point  D 
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Figure  74.  Coefficients  of  Sound  Speed  of  Point  D 
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Table  20.  EOF  Analysis  at  Point  D 


Mode  Number 

EOF  Mode 

PC 

Sound  Speed 

Depth 

Value 

Time 

Value 

1 

20 

Positive 

1960-1975 

1990-1992 

2000-2014 

Positive 

Increase 

80 

Positive 

Increase 

300 

Positive 

Increase 

800 

Zero 

No  change 

1 

20 

Positive 

1980-1990 

1992-2000 

Negative 

Decrease 

80 

Positive 

Decrease 

300 

Positive 

Decrease 

800 

Zero 

No  change 

2 

20 

Negative 

1960-1988 

2002-2014 

Positive 

Decrease 

80 

Negative 

Decrease 

300 

Positive 

Increase 

800 

Zero 

No  change 

2 

20 

Negative 

1988-2002 

Negative 

Increase 

80 

Negative 

Increase 

300 

Positive 

Decrease 

800 

Zero 

No  change 

3 

20 

Negative 

1960-1973 

1988-1998 

Positive 

Decrease 

80 

Positive 

Increase 

300 

Negative 

Decrease 

800 

Zero 

No  change 

3 

20 

Negative 

1998-2014 

Negative 

Increase 

80 

Positive 

Decrease 

300 

Negative 

Increase 

800 

Zero 

No  change 

4 

20 

Positive 

1960-1975 

1990-2012 

Positive 

Increase 

80 

Negative 

Decrease 

300 

Negative 

Decrease 

800 

Zero 

No  change 

4 

20 

Positive 

1975-1990 

Negative 

Decrease 

80 

Negative 

Increase 

300 

Negative 

Increase 

800 

Zero 

No  change 

At  Point  D  (Figures  73  and  74,  and  Table  20),  there  was  a  very  high  variability 
from  1990  to  2000  and  from  2000  to  2014.  The  sound  speed  increased  in  the  upper  layers 
from  1990  to  1992  and  from  2000  to  2014.  Conversely,  it  decreased  from  1992  to  2000  in 
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the  upper  layers.  Sound  speed  did  not  change  below  a  depth  of  800  meters  over  the  entire 
time  period. 

EOFs  of  Sound  Speed 


EOF-1  EOF-2  EOF -3  EOF-4 


-10  -5  0  -10  0  10  -5  0  5  -2  0  2 

Figure  75.  The  First  Four  EOFs  of  Point  E 


EOF  Coffecients  of  Sound  Speed 


1960  1970  I960  1990  2000  2010 

Year 

Figure  76.  Coefficients  of  Sound  Speed  of  Point  E 
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Table  21.  EOF  Analysis  at  Point  E 


Mode  Number 

EOF  Mode 

PC 

Sound  Speed 

Depth 

Value 

Time 

Value 

1 

5 

Negative 

1960-1970 

1970-1980 

1980-1990 

Positive 

Decrease 

25 

Negative 

Decrease 

40 

Negative 

Decrease 

70 

Negative 

Decrease 

1 

5 

Negative 

1990-2002 

2012-2014 

Negative 

Increase 

25 

Negative 

Increase 

40 

Negative 

Increase 

70 

Negative 

Increase 

2 

5 

Positive 

1992-2014 

Positive 

Increase 

25 

Positive 

Increase 

40 

Negative 

Decrease 

70 

Negative 

Decrease 

2 

5 

Positive 

1988-1992 

Negative 

Decrease 

25 

Positive 

Decrease 

40 

Negative 

Increase 

70 

Negative 

Increase 

3 

5 

Negative 

1973-2009 

Positive 

Decrease 

25 

Positive 

Increase 

40 

Positive 

Increase 

70 

Negative 

Decrease 

3 

5 

Negative 

1960-1972 

2010-2014 

Negative 

Increase 

25 

Positive 

Decrease 

40 

Positive 

Decrease 

70 

Negative 

Increase 

4 

5 

Negative 

1970-1980 

1990-2012 

Positive 

Increase 

25 

Positive 

Decrease 

40 

Negative 

Decrease 

70 

Positive 

Increase 

4 

5 

Negative 

1960-1970 

1980-1990 

1999-2014 

Negative 

Increase 

25 

Positive 

Decrease 

40 

Negative 

Increase 

70 

Positive 

Decrease 
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We  found  that  the  sound  speed  decreased  from  1960  to  1990  at  all  depths  at  point 
E  (Figures  75  and  76,  and  Table  21).  Conversely,  it  increased  from  1990  to  2002  and 
from  2012  to  2014.  There  was  a  very  high  variability  at  Point  E,  especially  from  1960  to 
1970  and  from  2000  to  2014. 


C.  TRANSMISSION  LOSS 

BELLHOP  shows  2-dimensional  ray  tracing  plots  for  the  given  sound  speed 
profiles  and  bottom  sediments.  In  this  section,  we  plotted  TL  that  BELLHOP  creates  by 
months  between  1960  and  2014  in  order  to  see  TL  variations. 


Figures  77-81  are  time -range  plots  of  TL  for  a  sound  receiver  located  at  50  m.  X- 
axis  is  time  labeled  in  years  and  have  a  1  month  resolution;  y-axis  is  range  along  the 
west-east  path  in  meters.  TL  values  are  mapped  as  pseudo-color  images  using  the  “jet” 
colormap.  The  “warm,”  yellow-orange  part  of  the  colorbar  corresponds  to  lower  TL 
values  and  higher  Received  Level.  The  “cold,”  blue-green  part  of  the  colorbar 
corresponds  to  higher  TL  and  lower  Received  Level. 


Figure  77.  Time-range  TL  plot  for  Point  A 
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Range,  m  Range, 


Point  B 


1960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010 


Figure  78.  Time-range  TL  plot  for  Point  B 


I960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010 

Figure  79.  Time-range  TL  plot  for  Point  C 
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Transmission  Loss,  d8  Transmission  Loss,  dB 


Range,  m  Range, 


Point  D 


1960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010 

Figure  80.  Time-range  TL  plot  for  Point  D 


I960  1965  1970  1975  1980  1985  1990  1995  2000  2005  2010 

Figure  8 1 .  Time-range  TL  plot  for  Point  E. 
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Transmission  Loss,  dB  Transmission  Loss,  dB 


At  point  A  (Figure  77),  TL  exhibits  a  high  inter-annual  variability  mostly  in 
magnitude.  The  highest  TL  at  long  ranges  occurred  between  1982  and  1990.  Time 
periods  from  1970  to  1982  and  from  2010  to  2014  seemed  to  have  very  low  TLs  overall. 
Conversely,  the  time  span  from  1982  to  1990  at  close  ranges,  0-5  km,  showed  a  lower 
TL  than  from  1970  to  1982. 

Point  B  (Figure  78)  represented  shallow  water  in  the  Central  Mediterranean  Sea. 
At  this  location,  the  main  patterns  of  acoustic  propagation  are  bottom  bounce  and  surface 
duct.  Figure  79  does  not  show  as  much  inter-annual  variability  as  Point  A.  This  may  be  a 
result  of  the  shallow  water  and  interactions  with  bottom. 

At  point  C  (Figure  79),  TL  changed  quickly  between  1960  and  1995.  After  1995, 
the  TL  exhibited  the  same  values.  Thus,  we  assume  that  the  TL  after  1995  shows  a 
different  TL  distribution  than  before  1995.  In  addition,  the  TL  appeared  lower  between 
1995  and  2014.  Thus,  we  assume  that  detection  ranges  increased  after  1990.  Horizontal 
TL  lines  represent  convergence  zone  paths.  Therefore,  we  found  that  a  lower  TL  might 
occur  at  longer  ranges  due  to  the  presence  of  caustics  at  the  convergence  zone. 

Three  horizontal  TL  lines  in  Figure  80  represent  the  convergence  zones  at  the 
point  D  location.  They  occur  at  the  same  range  over  the  all  decades,  and  they  do  not  show 
high  inter- annual  variability.  Ranges  closer  than  38  km  show  high  inter-annual 
variability.  There  were  low  TLs  from  1960  to  1978,  from  1985  to  1992,  and  from  2007  to 
2014.  Conversely,  high  TLs  manifested  from  1978  to  1985  and  from  1992  to  2007.  The 
lowest  TL  occurred  between  2007  and  2014. 

Three  horizontal  TL  lines  occur  at  the  same  range  over  all  decades  in  Figure  81 
for  the  point  E  location.  They  do  not  show  a  high  variability.  Ranges  closer  than  38  km 
show  inter-annual  variability;  however,  it  takes  place  periodically.  Thus,  TLs  could  not 
be  categorized  into  decades.  Nevertheless,  the  lowest  TL  occurred  between  1967 
and  1973. 
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D.  TRANSMISSION  LOSS  DUE  TO  SHORT-TERM  INTER-ANNUAL 
VARIABILITY 

August  and  January  have  been  selected  to  represent  the  summer  and  winter 
seasons.  We  selected  time  periods  according  to  years  with  very  different  TL  distributions 
to  show  the  effect  of  short-term  inter-annual  variability. 


a.  Point  A 


In  the  previous  section,  we  saw  that  Point  A  had  a  lower  TL  in  the  period  from 
1970  to  1982  and  from  2010  to  2014.  In  order  to  see  the  TL  variations  between  these  time 
periods,  TL  was  modeled  using  BELLHOP  with  sound  speed  profiles  averages  over  the 
selected  intervals.  Results  are  plotted  in  Figure  82. 
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Figure  82.  TL  for  1970-1979  (upper  panels),  1980-1990  (middle  panels) 
and  2000-2014  (lower  panels)  for  Point  A  in  January 
(left  column)  and  August  (right  column) 
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Figure  82  (cont’d)  TL  for  1970-1979  (upper  panels),  1980-1990  (middle  panels) 
and  2000-2014  (lower  panels)  for  Point  A  in  January 
(left  column)  and  August  (right  column) 


TL  variations  in  January  were  distributed  using  different  limits.  The  1970-1979  time 
period  (upper  left)  has  a  60  dB  TL  at  close  distances;  however,  it  has  a  TL  of  about  70- 
85  dB  at  approximately  3  km  and  45  km.  The  1980-1990  time  period  (middle  left)  shows 
a  TL  of  75-85  dB  between  the  same  ranges.  Specifically,  it  exhibits  a  lower  TL  at  the 
surface  duct  and  in  the  upper  ray  limit  of  the  convergence  zone.  As  compared  to  the 
2000-2014  time  period  (lower  left)  it  has  the  lower  TL  for  close  ranges,  but  it  has  a 
higher  TL  at  long  ranges.  These  results  correspond  with  what  we  found  in  the  previous 
section. 

TL  distributions  in  August  1970-1979  (upper  right)  and  1980-1990  (middle  right) 
show  the  same  TL  at  closer  ranges;  however,  the  1970-1979  time  period  shows  a  lower 
TL  at  long  ranges.  The  plot  for  August  2000-2014  (lower  right)  shows  the  same  TL  at 
closer  ranges  as  in  1980-1990  but  a  lower  TL  at  long  distances. 

The  researchers  found  that  1970-1979  and  2000-2014  have  higher  TLs  than 
1980-1990  for  closer  ranges  in  January;  however,  they  have  the  same  TL  in  August.  For 
longer  ranges,  1980-1990  has  the  highest  TL.  Also,  we  found  that  January  shows 
significantly  more  variations  than  August. 
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b.  Point  B 

From  the  EOF  analysis,  we  saw  that  TL  does  not  change  very  much  inter- 
annually.  We  found  that  2007-2008  and  1997-1998  exhibited  the  greatest  variability. 
Sound  speed  throughout  all  depths  decreased  during  these  periods.  Figure  83  compares 
TL  for  1997-1998  (upper  panels)  and  2007-2008  (middle  panels)  with  the  SMG-WOD 
multi-year  monthly  averages  calculated  over  the  total  time  period  of  1960-2014  for 
January  (left  columns)  and  August  (right  columns). 

In  January,  the  surface  duct  weakens  and  bottom  bounce  strengthens  both 
between  1997  and  1998  (upper  left)  and  between  2007  and  2008  (middle  left),  as 
compared  to  the  multi-year  average  (lower  left).  In  August,  the  bottom  bounce  weakens 
and  changes  into  a  surface  sound  channel  between  1997  and  1998  (upper  right).  August 
1997-1998  also  provides  lower  TL  at  longer  ranges.  In  August  2007-2008  (middle  right) 
the  bottom  bounce  does  not  change  much,  a  higher  TL  occurred  in  all  ranges,  and  there 
was  greater  propagation  at  the  surface. 
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Figure  83.  TL  for  1997-1998  (upper  panels),  2007-2008  (middle  panels) 
and  multi-year  1960-2014  monthly  averages  (lower  panels)  for 
Point  B  in  January  (left  column)  and  August  (right  column) 
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Figure  83  (cont’d.).  TL  for  1997-1998  (upper  panels),  2007-2008  (middle  panels) 
and  multi-year  1960-2014  monthly  averages  (lower  panels)  for  Point  B  in  January 
(left  column)  and  August  (right  column) 


These  results  confirmed  that  there  was  a  high  variability  from  1997  to  1998  and 
from  2007  to  2008.  Also,  we  confirmed  that  January  shows  a  lot  more  variations 
than  August. 


c.  Point  C 

We  found  that  there  is  very  high  variability  from  1980  to  1995  and  from  2000  to 
2014  based  on  EOF  analysis.  Also,  the  TL  after  1995  exhibited  a  different  distribution 
than  before  1995.  For  these  reasons,  we  plotted  the  TL  for  1980-1995  and  2000-2014 
(Figure  84)  in  January  (left  column)  and  August  (right  column). 
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Figure  84.  TL  for  1980-1995  (upper  panels)  and  2000-2014  (lower  panels) 
for  Point  C  in  January  (left  column)  and  August  (right  column) 


There  are  some  important  variations  of  TL  between  these  years.  The  1980-1995 
time  period  in  January  (upper  left)  only  shows  convergence  propagation  and  the  upper 
limiting  ray  of  the  convergence  zone  that  extends  to  a  depth  of  800-1000  meters. 
However,  the  right  figure  shows  a  surface  duct  and  convergence  zone  propagation 
together.  Also,  in  2000-2014  time  period  in  January  (lower  left),  the  upper  limit  of  the 
convergence  zone  extends  to  a  depth  of  900-1100  meters. 

In  August,  both  the  time  periods  show  convergence  zone  propagation.  However, 
there  is  a  slight  difference  between  the  upper  limiting  ray  of  the  convergence  zone.  The 
upper  limiting  ray  extends  to  a  depth  of  1500  meters  in  1980-1995  and  a  depth  of  1600 
meters  in  2000-2014.  Also,  January  showed  a  lot  more  variations  than  August.  This 
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anomaly  must  have  been  the  result  of  mixing  layer  depths  and  sound-speed  profiles 
because  the  Mediterranean  Sea  has  very  strong,  severe  winds  in  winter,  but  not  in 
summer.  These  weather  patterns  may  create  more  inter-annual  variations  at  the  surface, 
which  affect  sound-speed  profiles.  Although  we  saw  some  inter-annual  variations  in 
August,  these  actually  did  not  affect  sound-speed  profiles  very  much.  The  variations  may 
be  a  result  of  very  warm  water  at  the  surface.  Despite  the  inter-annual  variations,  August 
exhibited  the  highest  sound  speed  at  the  surface  without  an  upward  bend. 

d.  Point  D 

From  the  EOF  analysis,  the  researchers  found  that  the  years  from  1990  to  2000 
and  from  2000  to  2014  exhibited  high  variability.  For  this  reason,  we  plotted  the  TL  for 
1990-2000  and  2000-2014  (Figure  85). 

The  sound  propagation  profiles  in  January  were  almost  the  same  (Figure  85,  left 
column).  Although  there  was  high  variability  between  these  years,  it  did  not  affect  the 
sound  propagation  type  and  TL  very  much.  There  was  only  a  slight  difference  in  the 
shadow  zones  at  the  surface.  Thus,  we  found  there  was  a  slight  inter-annual  variability  at 
Point  D  in  January.  In  August,  both  the  time  intervals  demonstrate  almost  the  same  sound 
propagation  pattern  and  TL  (right  column).  The  1990-2000  propagation  shows  a  slightly 
lower  TL  than  in  2000-2014  at  long  ranges.  Therefore,  we  expected  a  slight  inter-annual 
variability  at  Point  D  in  August. 
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Figure  85.  TL  for  1990-2000  (upper  panels)  and  2000-2014  (lower  panels) 
for  Point  D  in  January  (left  column)  and  August  (right  column) 


e.  Point  E 

From  EOF  analysis,  we  found  that  Point  E  demonstrated  very  high  variability 
from  1960  to  1980  and  from  2000  to  2014.  We  also  found  that  it  showed  a  lower  TL 
between  1967  and  1973  from  TL  analysis.  In  order  to  see  inter-annual  variability  between 
these  spans,  we  plotted  Figure  86. 
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Figure  86.  TL  for  1960-1970  (upper  panels)  and  2000-2014  (lower  panels) 
for  Point  E  in  January  (left  column)  and  August  (right  column) 


In  January,  there  was  high  variability  between  1960-1970  and  2000-2014  (Figure 
86,  left  column).  The  earlier  interval  shows  a  lower  TL  at  the  surface,  propagated  over  a 
long  range  and  concentrated  between  the  surface  and  a  depth  of  300  meters.  The  later 
interval  demonstrates  a  sound  channel  and  surface  propagation  together  with  a  higher  TL 
at  long  ranges  and  a  lower  TL  at  short  ranges.  Therefore,  we  confirmed  that  Point  E 
showed  high  variability  in  January.  This  variability  must  have  been  the  result  of  high-rate 
sound-speed  changes  at  all  depths. 

In  August,  the  sound  propagations  are  almost  identical  in  1960-1970  and  2000- 
2014  (Figure  86,  right  column).  Both  the  intervals  demonstrate  surface  propagation; 
however,  the  earlier  one  shows  propagation  to  deeper  layers  than  the  later  one.  There  is 
also  a  difference  in  TLs  between  the  figures.  The  left  figure  shows  a  higher  TL  at  short 
ranges  and  a  lower  TL  at  long  ranges.  The  researchers  found  that  Point  E  had  a  higher 
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variability  in  January  than  in  August.  Moreover,  1960-1970  shows  a  lower  TL  at  short 
ranges  and  a  higher  TL  at  long  ranges  in  both  months.  We  also  examined  why  shallow 
water  shows  more  inter-annual  variations  than  deeper  waters.  The  findings  indicate  that 
the  mixed  layer,  especially  in  winter,  may  show  higher  variability  than  deep  layers. 


120 


VII.  CONCLUSIONS 


The  primary  area  of  research  was  inter-annual  variability  of  sound  propagation  in 
the  Mediterranean  Sea  using  synoptic  monthly  gridded  data  and  environmental 
characteristics  from  the  Generalized  Digital  Environmental  Model  (GDEM).  The  latter 
has  been  widely  used  in  to  model  acoustic  propagation  for  Navy  needs,  while  the  former 
is  relatively  new  and  has  not  been  applied  previously  for  acoustic  modeling  in  the 
Mediterranean  Sea. 

We  confirmed  that  both  the  climatological  and  multi-year  mean  synoptic 
databases  correctly  reproduce  the  high  stratification  of  the  Mediterranean  Sea  and  its 
characteristic  water  masses.  In  particular,  we  found  that  the  Levantine  Surface  Water 
(LSW),  which  forms  in  August,  has  the  highest  spiciness  (9. 1  m3 /kg)  among  all  water 
masses  in  the  Mediterranean  Sea.  The  Eastern  Mediterranean  Sea  shows  a  greater  scale  of 
variability  of  spiciness  than  the  Western  Mediterranean  Sea  in  August,  but  not  in  January. 
We  also  found  that  spiciness  decreases  in  summer  and  increases  in  winter.  Both  basins 
showed  the  same  scale  of  variation  in  January,  but  not  in  August.  In  addition,  there  are 
high  spice  values  around  the  Crete  passage  where  the  Levantine  Intermediate  Water 
(LIW)  forms  in  winter  in  response  to  strong  northerly  winds.  High  salinity  and 
temperature,  or  high  spiciness,  is  the  signature  of  the  LIW. 

Five  locations  have  been  selected  in  the  Western,  Central  and  Eastern 
Mediterranean  basins  to  reflect  spatial  differences  in  seasonal  and  inter-annual 
variability,  as  well  as  account  for  different  geoacoustic  characteristics  and  bathymetry.  A 
composite  analysis  have  been  applied  to  sound  speed  profiles  at  each  locations.  This 
allowed  to  separate  the  seasonal  and  inter-annual  variability  in  the  Mediterranean  sea. 

The  temperature  profiles  at  all  the  locations  exhibit  strong  seasonality  associated 
with  surface  heat  fluxes  and  variations  in  the  vertical  structure  of  water  column. 
However,  this  seasonality  was  strongly  modulated  by  inter-annual  and  decadal 
variability.  For  salinity  the  inter-annual  variability  is  much  stronger  than  seasonal,  most 
likely  associated  with  inter-annual  and  decadal  variability  of  freshwater  fluxes  and  water 
mass  formation.  For  sound  speed,  variability  shows  that  the  upper  layer  sound  speed  is 
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determined  by  temperature  as  expected,  while  the  lower  level  is  determined  by  the 
pressure  and  shows  neither  seasonal  nor  inter-annual  variability.  It  is  however  expected 
that  strong  inter-annual  variability  of  salinity  will  affect  the  inter-annual  variability  of 
sound  speed.  Although  temperature  and  salinity  decrease  or  increase  in  some  decades, 
overall,  the  Mediterranean  is  becoming  warmer  and  saltier  from  1960  to  2014,  especially 
from  2000  to  2014. 

The  SMG-WOD  multi-year  monthly  averages  agreed  well  with  and  GDEM 
monthly  mean  sound-speed  profdes.  There  was  a  slight  difference  between  them  in  the 
upper  layers,  especially  in  the  thermoclines.  There  are  two  possible  reasons.  First,  the 
vertical  resolution  is  different  in  these  databases.  In  the  upper  part  of  the  Mediterranean 
Sea,  GDEM  has  78  vertical  layers  of  two-meter  vertical  resolution;  however,  SMG-WOD 
has  57  vertical  layers  with  a  5  m  step.  Second,  GDEM  means  are  calculated  for  1945— 
2014;  however,  the  SMG-WOD  covers  data  between  1960  and  2014.  This  15-years’ 
difference  can  result  in  biased  mean  values  due  to  significant  inter-annual  and  decadal 
variability  in  the  Mediterranean  Sea 

Empirical  orthogonal  functions  (EOF)  analysis  was  conducted  to  identify  the 
inter-annual  variability  of  the  sound  speed  profiles.  The  number  of  modes  necessary  to 
capture  90%  of  total  variance  depends  on  the  region.  Point  E  is  very  shallow  and  needs 
only  three  modes  to  reach  90  percent.  In  the  Eastern  Mediterranean  Sea,  Point  A  and 
Point  C  demand  10-11  modes  to  reach  a  90  percent  variance;  however,  Point  D  in  the 
Central  Mediterranean  Sea  reaches  90  percent  with  only  six  modes.  Point  B,  which  is 
shallower,  reaches  90  percent  with  10  modes.  Thus,  we  cannot  generalize  which  region 
demands  more  modes.  We  can  only  confirm  that  Points  A  and  C  in  the  Eastern 
Mediterranean  need  more  modes  than  Point  D  at  the  same  depth  in  the  Central 
Mediterranean. 

The  maximum  sound-speed  variability  was  observed  in  different  decades  for 

different  points.  Point  A  (a  deep  location  in  the  Eastern  basin)  had  great  variability  from 

1980  to  1990  and  from  2000  to  2014.  The  largest  variability  was  observed  in  2009-2010. 

Point  B  (central  shallow  location)  exhibited  the  greatest  variability  in  two  two-year 

periods:  2007-2008  and  1997-1998.  Point  C  (another  deep  location  in  the  Eastern 
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Mediterranean)  had  high  variability  in  1980-1995  and  2000-2014.  Point  D,  which  is 
located  in  deep  water  of  the  central  part,  shows  great  variability  in  1990-2014.  Point  E  is 
in  very  shallow  water  of  the  Western  Mediterranean  Sea,  and  has  great  variability  in  all 
years.  The  maximum  variations  took  place  between  2000  and  2014. 

The  inter-annual  variability  detected  by  the  EOF  analysis  may  result  in  different 
shapes  of  the  sound  speed  profiles  in  different  years  for  the  same  seasons.  To  assess  how 
this  variability  affects  acoustic  propagation  in  the  Mediterranean  sea,  a  ray  trace  model 
Bellhop  was  used.  The  sound  source  was  located  at  40  m  and  has  frequency  of  3500  Hz. 
The  propagation  was  modeled  in  west-east  direction  up  to  a  maximum  range  of  70 
nautical  miles.  TL  profiles  modeled  with  SMG-WOD  multi-year  monthly  average  are 
similar  to  those  produced  by  GDEM  monthly  mean  sound  speed  profiles.  Nevertheless, 
there  are  some  differences.  For  example,  the  lowest  limiting  ray  of  the  convergence  zone 
extends  deeper  in  the  SMG-WOD  than  the  GDEM  in  very  deep  layers. 

TL  was  calculated  for  the  five  selected  locations  with  the  acoustic  environment 
modeled  by  averaging  sound  speed  profiles  over  different  time  intervals.  These  time 
intervals  were  chosen  based  on  the  results  of  the  EOF  analysis.  The  acoustic  modeling 
was  conducted  for  January  and  August  since  it  was  determined  that  the  maximum 
variance  of  sound  speed  in  the  Eastern  Mediterranean  Sea  occurs  from  January  to  March; 
however,  for  the  Central  and  Western  Mediterranean  Sea,  the  maximum  variance  occurs 
in  August. 

We  analyzed  TLs  for  each  decade  and  found  that  the  TL  did  not  show  high  inter¬ 
annual  variability  over  convergence  zone  paths.  We  also  found  that  TL  showed  higher 
inter-annual  variability  at  short,  up  to  5  km  ranges,  but  not  for  long  ranges.  January 
showed  more  inter-annual  changes  than  August.  We  detected  different  types  of  sound 
propagation  in  January,  but  not  in  August.  A  more  significant  inter-seasonal  variability 
was  observed  in  the  upper  layer  of  the  water  column.  We  hypothesized  that  it  may  be  due 
greater  variability  of  the  mixed  layer  as  compared  to  deep  layers,  especially  in  the  winter. 

The  observed  inter-annual  variability  in  the  acoustic  propagation  can  be  important 
for  tactical  naval  applications.  The  convergence  zones  occur  in  the  deep  layer  of  the 
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Mediterranean  Sea  in  all  seasons.  Thus,  it  is  always  possible  to  detect  submarines  from 
long  distances  at  the  surface.  Furthermore,  a  surface  duct  emerges  in  winter  and 
propagates  sound  over  long  distances.  Therefore,  we  can  say  that  the  winter  season 
creates  a  disadvantageous  environment  for  submarines  sailing  at  periscope  depths. 
Although  both  the  datasets  exhibit  similar  TL  profiles,  there  are  some  significant 
differences  that  can  affect  naval  operations.  For  example,  the  SMG-WOD  showed  an 
uninterrupted  surface  duct  at  Point  A  in  January;  however,  the  GDEM  exhibited  an 
interrupted  surface  duct,  which  suggests  that  submarines  cannot  be  detected  from  the 
surface  at  some  specific  ranges.  Additionally,  there  are  few  differences  in  the  TL  values 
of  the  shadow  zones.  Shadow  zones  are  very  important  regions  for  submarines  to  hide. 
Submarines  can  select  shadow  zones  where  ranges  are  poor  to  avoid  detection. 
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APPENDIX  A.  BOTTOM  SEDIMENTS 


Table  22.  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

PT 

Sand 

102 

1125 

25 

9 

PT 

Silt 

105 

1105 

52 

18 

PT 

Silty  Clay 

107 

1107 

50 

7 

PT 

Clay 

108 

1108 

8 

23 

PT 

Gravel 

112 

1112 

62 

3 

PT 

Sandy  Mud 

196 

1104 

25 

96 

PT 

Mud 

197 

1107 

7 

PC 

Sand 

202 

1225 

25 

9 

PC 

Silty  Sand 

203 

1203 

24 

11 

PC 

Sandy  Silt 

204 

1204 

53 

16 

PC 

Silt 

205 

1205 

52 

18 

PC 

Clayey  Silt 

206 

1206 

19 

PC 

Silty  Clay 

207 

1207 

50 

22 

PC 

Clay 

208 

1208 

8 

23 

PC 

Sand  -  Silt  -  Clay 

209 

1209 

52 

17 

PC 

Marl 

211 

1208 

8 

23 

PC 

Gravel 

212 

1212 

62 

3 

PC 

Gravelly  Sand 

220 

1220 

21 

7 

PC 

Coarse  Sand 

222 

1222 

22 

7 

PC 

Medium  Sand 

223 

1223 

23 

9 

PC 

242 

1242 

25 

20 

PC 

Ooze 

310 

1208 

8 

23 

PS 

Silt 

405 

1105 

52 

18 

PS 

Silty  Clay 

7 

PS 

Clay 

408 

1108 

8 

23 

PS 

Ooze 

510 

1108 

8 

23 

PS 

Mud 

511 

1105 

52 

17 

p 

Clayey  Silt 

706 

1106 

50 

19 

p 

Silty  Clay 

1107 

50 

22 

p 

Clay 

708 

1108 

8 

23 

V 

Sand 

802 

1102 

23 

9 

V 

Silt 

805 

1105 

52 

18 

V 

Silty  Clay 

807 

50 

7 

o 

NO  DATA 

888 

888 

888 

888 

L 

LAND 

999 

999 

999 

999 

r  T 

Rock 

1101 

1101 

1 

2 

T 

Sand 

1102 

23 

9 

T 

Silty  Sand 

1103 

25 

11 

T 

Sandy  Silt 

25 

16 

T  ' 

Silt 

1105 

1105 

52 

18 

125 


Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

T 

Clayey  Silt 

1106 

1106 

19 

T 

Silty  Clay 

1107 

22 

T 

Clay 

1108 

8 

23 

T 

Sand  -  Silt  -  Clay 

1109 

24 

17 

T 

Gravel 

1112 

1112 

62 

3 

T 

Sandy  Gravel 

1113 

1113 

63 

4 

T 

Silty  Gravel 

1114 

1114 

63 

10 

T 

Muddy  Sandy  Gravel 

1115 

1114 

63 

6 

T 

Clayey  Gravel 

1116 

1116 

63 

T 

Muddy  Gravel 

1117 

1116 

63 

T 

Gravelly  Muddy  Sand 

1118 

1120 

63 

8 

T 

Gravel  -  Silty  Sand 

1119 

1113 

63 

8 

T 

Gravelly  Sand 

21 

7 

T 

Very  Coarse  Sand 

1121 

1121 

21 

5 

T 

Coarse  Sand 

1122 

1122 

22 

7 

T 

Medium  Sand 

1123 

1123 

23 

9 

T 

Fine  Sand 

1124 

1124 

24 

11 

T 

Very  Fine  Sand 

1125 

1125 

25 

13 

T 

Clayey  Sand 

1126 

1126 

25 

14 

T 

Gravel  -  Shell 

1128 

1112 

62 

3 

T 

Gravelly  Silt 

1130 

1130 

21 

16 

T 

Gravelly  Silt  -  Shell 

1131 

21 

16 

T 

Gravelly  Sandy  Silt 

1132 

1130 

21 

16 

T 

Gravelly  Mud 

1133 

21 

16 

T 

Gravel  -  Sand  -  Mud 

1134 

1113 

63 

6 

T 

Rock  -  Sand  -  Mud 

1135 

1161 

61 

3 

T 

Rock  -  Gravel  -  Mud 

1136 

1161 

61 

3 

T 

Rock  -  Gravel  -  Sand 

1137 

1161 

61 

3 

T 

Rock  -  Gravel  -  Sand  -  Mud 

1138 

1161 

61 

3 

T 

Gravelly  Clay 

1140 

1140 

21 

21 

T 

Sand  -  Clay  -  Shell 

1141 

23 

9 

T 

Sandy  Clay 

1142 

1142 

25 

20 

T 

Coarse  Sand  -  Shell 

1146 

1113 

63 

4 

T 

Very  Fine  Silt 

1150 

50 

21 

T 

Fine  Silt 

1151 

1151 

51 

19 

T 

Medium  Silt 

1152 

1152 

52 

17 

T 

Coarse  Silt 

1153 

1153 

53 

15 

T 

Rough  Rock 

1154 

1101 

1 

1 

T 

Mud  over  Rock 

1155 

1101 

1 

2 

T 

Silty  Clay  -  Shell 

1156 

1130 

21 

16 

T 

Boulders 

1160 

1160 

60 

2 

T 

Cobbles  (Stones)  -  Shell 

1161 

1161 

61 

3 

T 

Pebbles  -  Shell 

1162 

1162 

62 

3 

T 

Granules 

1163 

1163 

63 

3 

T 

Sand  -  Silt  -  Clay  -  Shell 

1164 

1112 

62 

16 

T 

Gravel  -  Sand  -  Shell 

1165 

1161 

61 

3 

T 

Shell 

1166 

1162 

62 

3 

126 


Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

c 

Muddy  Gravel 

1217 

1216 

63 

10 

c 

Gravelly  Muddy  Sand 

1218 

1220 

63 

8 

c 

Gravel  -  Silty  Sand 

1219 

1213 

63 

8 

c 

Gravelly  Sand 

1220 

21 

7 

c 

Very  Coarse  Sand 

1221 

1221 

21 

5 

c 

Coarse  Sand 

1222 

1222 

22 

7 

c 

Medium  Sand 

1223 

1223 

23 

9 

c 

Fine  Sand 

1224 

1224 

24 

11 

c 

Very  Fine  Sand 

1225 

1225 

25 

13 

c 

Clayey  Sand 

1226 

1226 

25 

14 

c 

Oolite 

1227 

1223 

23 

9 

c 

Gravel  -  Shell 

1228 

1212 

62 

3 

c 

Gravelly  Silt 

1230 

1230 

21 

16 

c 

Gravelly  Silt  -  Shell 

1231 

1230 

21 

16 

c 

Gravelly  Sandy  Silt 

1232 

1230 

21 

16 

c 

Gravelly  Mud 

1233 

1230 

21 

16 

c 

Gravel  -  Sand  -  Mud 

1234 

1213 

63 

6 

c 

Rock  -  Sand  -  Mud 

1235 

1261 

61 

3 

c 

Rock  -  Gravel  -  Mud 

1236 

1261 

61 

3 

c 

Rock  -  Gravel  -  Sand 

1237 

1261 

61 

3 

c 

Rock  -  Gravel  -  Sand  -  Mud 

1238 

1261 

61 

3 

c 

Gravelly  Clay 

1240 

21 

21 

c 

Sand  -  Clay  -  Shell 

1241 

23 

9 

c 

inM  53  Ul«ll 

1242 

1242 

25 

20 

c 

Coral  Debris  -  Sand 

1243 

1212 

62 

3 

c 

Coral  Debris  -  Sand  -  Shell 

1244 

1212 

62 

3 

c 

Coral  Debris  -  Shell 

1245 

1212 

62 

3 

c 

Coarse  Sand  -  Shell 

1246 

1213 

63 

4 

c 

Coral  Debris  -  Sand  -  Mud 

1247 

1212 

62 

6 

c 

Coral  Debris  -  Mud  -  Shell 

1248 

1212 

62 

6 

c 

Coral  Debris  -  Mud 

1249 

1212 

62 

6 

c 

Very  Fine  Silt 

1250 

1250 

50 

21 

c 

Fine  Silt 

1251 

1251 

51 

19 

c 

Medium  Silt 

1252 

1252 

52 

17 

c 

Coarse  Silt 

1253 

1253 

53 

15 

c 

Rough  Rock 

1254 

1 

1 

c 

Mud  over  Rock 

1255 

1201 

1 

2 

c 

Silty  Clay  -  Shell 

1256 

1230 

21 

16 

Coral  Debris  -  Sand  -  Mud  - 
Shell 

1257 

1212 

62 

3 

c 

Coral  Debris 

1258 

1212 

62 

3 

c 

Coral 

1259 

1201 

1 

1 

c 

Boulders 

1260 

1260 

60 

2 

c 

Cobbles  (Stones)  -  Shell 

1261 

1261 

61 

3 

c 

Pebbles  -  Shell 

1262 

1262 

62 

3 

c 

Granules 

1263 

1263 

63 

3 

c 

Sand  -  Silt  -  Clay  -  Shell 

1264 

1212 

62 

16 

127 
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Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


pe 

Enhanced  Categories 
Description 

Standard 

Reduced 

Gravel  -  Sand  -  Shell 

1265 

1261 

61 

3 

Shell 

1266 

1262 

62 

3 

Rock  -  Gravel  -  Sand  -  Shell 

1267 

1261 

61 

3 

Sand  -  Shell 

1268 

1213 

63 

4 

Rock  -  Gravel 

1270 

1201 

1 

2 

Rock  -  Coral 

1271 

1201 

1 

1 

Rock  -  Sand 

1272 

1201 

1 

2 

Rock  -  Mud 

1273 

1201 

1 

2 

Mud  -  Shell 

1274 

1214 

63 

10 

Gravel  -  Sand 

1275 

1213 

63 

4 

Gravel  -  Mud 

1276 

1214 

63 

10 

Clayey  Sand  -  Shell 

1277 

1220 

21 

7 

Soft  Mud 

1278 

1206 

50 

21 

Hard  Mud 

1279 

1206 

50 

19 

Silty  Sand  -  Shell 

1281 

1221 

21 

5 

Gravelly  Sand  -  Shell 

1282 

1213 

63 

4 

Medium  Sand  -  Shell 

1283 

1220 

21 

5 

Fine  Sand  -  Shell 

1284 

1220 

21 

7 

Sandy  Gravel  -  Shell 

1285 

1213 

63 

4 

Clayey  Silt  -  Shell 

1286 

1230 

21 

16 

Silt  -  Shell 

1287 

1230 

21 

16 

Silty  Gravel  -  Shell 

1288 

1214 

63 

10 

Sandy  Silt  -  Shell 

1289 

1213 

21 

8 

Muddy  Tidal  Flats 

1253 

53 

15 

Sandy  Tidal  Flats 

1291 

1223 

23 

9 

Sandy  Muddy  Tidal  Flats 

1292 

1209 

24 

17 

Sand  Dune 

1293 

1202 

23 

9 

Mud  -  Sand 

1294 

1209 

24 

17 

Muddy  Sand 

1295 

1209 

24 

12 

Sandy  Mud 

1296 

1204 

25 

18 

Mud 

1297 

1205 

52 

18 

Clay  -  Shell 

1298 

21 

10 

Stiff  Mud 

1299 

50 

19 

Silt 

1405 

52 

18 

Clay 

1108 

8 

23 

Ooze 

1108 

8 

23 

Mud 


Sand 


Silt 


Gravel 


Sandy  Gravel 


Gravelly  Sand 


Rough  Rock 


Boulders 


Rock  -  Gravel 


Rock  -  Sand 


Gravel  -  Sand 


Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

uc 

Silt 

2205 

1205 

52 

18 

uc 

Clayey  Silt 

2206 

1206 

50 

19 

uc 

Silty  Clay 

2207 

1207 

50 

22 

uc 

Clay 

2208 

1208 

8 

23 

uc 

Sand  -  Silt  -  Clay 

2209 

1209 

24 

17 

uc 

Marl 

2211 

1208 

8 

23 

US 

Silt 

2405 

1105 

52 

18 

US 

Clay 

2408 

1108 

8 

23 

uc 

Ooze 

2210 

1208 

8 

23 

u 

Silty  Sand 

2603 

25 

11 

u 

Sandy  Silt 

2604 

25 

16 

u 

Silt 

2605 

52 

18 

u 

Clayey  Silt 

2606 

1106 

50 

19 

u 

Silty  Clay 

2607 

1107 

50 

22 

u 

Clay 

2608 

1108 

8 

23 

u 

Sand  -  Silt  -  Clay 

2609 

1109 

24 

17 

u 

Fine  Sand 

2624 

1124 

24 

11 

u 

Gravel  -  Sand 

2675 

1113 

63 

4 

u 

Gravelly  Sand  -  Shell 

2682 

1113 

63 

5 

HC 

Rock 

3201 

1201 

1 

2 

HC 

Sand 

3202 

1202 

23 

9 

HC 

Silty  Sand 

3203 

1203 

25 

11 

HC 

Sandy  Silt 

3204 

1204 

25 

16 

HC 

Silt 

3205 

1205 

52 

18 

HC 

Clayey  Silt 

3206 

1206 

50 

19 

HC 

Silty  Clay 

3207 

1207 

50 

22 

HC 

Clay  (Marl) 

3208 

1208 

8 

23 

HC 

Sand  -  Silt  -  Clay 

3209 

1209 

24 

17 

HC 

Ooze 

3210 

1208 

8 

23 

HC 

Gravel  (Shell  Detritus) 

3212 

1212 

62 

3 

HC 

Sandy  Gravel 

3213 

1213 

63 

4 

HC 

Silty  Gravel 

3214 

1214 

63 

10 

HC 

Muddy  Sandy  Gravel 

3215 

1214 

63 

6 

HC 

Clayey  Gravel 

3216 

1216 

63 

10 

HC 

Muddy  Gravel 

3217 

1216 

63 

10 

HC 

Gravelly  Muddy  Sand 

3218 

1220 

63 

8 

HC 

Gravel  -  Silty  Sand 

3219 

1213 

63 

8 

HC 

Gravelly  Sand 

3220 

1220 

21 

7 

HC 

Very  Coarse  Sand 

3221 

1221 

21 

5 

HC 

Coarse  Sand 

3222 

1222 

22 

7 

HC 

Medium  Sand 

3223 

1223 

23 

9 

HC 

Fine  Sand 

3224 

1224 

24 

11 

HC 

Very  Fine  Sand 

3225 

1225 

25 

13 

HC 

Clayey  Sand 

3226 

1226 

25 

14 

HC 

Gravel  -  Shell 

3228 

1212 

62 

3 

HC 

Gravelly  Silt 

3230 

1230 

21 

16 

HC 

Gravelly  Silt  -  Shell 

3231 

1230 

21 

16 
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Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

HC 

Gravelly  Sandy  Silt 

3232 

1230 

21 

16 

HC 

Gravelly  Mud 

3233 

1230 

21 

16 

HC 

Gravel  -  Sand  -  Mud 

3234 

1213 

63 

6 

HC 

Rock  -  Sand  -  Mud 

3235 

1261 

61 

3 

HC 

Rock  -  Gravel  -  Mud 

3236 

1261 

61 

3 

HC 

Rock  -  Gravel  -  Sand 

3237 

1261 

61 

3 

HC 

Rock  -  Gravel  -  Sand  -  Mud 

3238 

1261 

61 

3 

HC 

1240 

21 

21 

HC 

Sand  -  Clay  -  Shell 

3241 

1202 

23 

9 

HC 

Sandy  Clay  (Sandy  Marl) 

3242 

1242 

25 

HC 

Coral  Debris  -  Sand 

3243 

1212 

62 

3 

HC 

Coral  Debris  -  Sand  -  Shell 

3244 

1212 

62 

3 

HC 

Coral  Debris  -  Shell 

3245 

1212 

62 

3 

HC 

Coarse  Sand  -  Shell 

3246 

1213 

63 

4 

HC 

Coral  Debris  -  Sand  -  Mud 

3247 

1212 

62 

6 

HC 

Coral  Debris  -  Mud  -  Shell 

3248 

1212 

62 

6 

HC 

Coral  Debris  -  Mud 

3249 

1212 

62 

6 

HC 

Very  Fine  Silt 

1250 

50 

21 

HC 

Fine  Silt 

3251 

1251 

51 

19 

HC 

Medium  Silt 

3252 

1252 

52 

17 

HC 

Coarse  Silt 

3253 

1253 

53 

15 

HC 

Rough  rock 

3254 

1201 

1 

1 

HC 

Mud  over  Rock 

3255 

1201 

1 

2 

HC 

Silty  Clay  -  Shell 

3256 

1230 

21 

16 

HC 

Coral  Debris  -  Sand  -  Mud  - 
Shell 

3257 

1212 

62 

3 

HC 

Coral  Debris 

3258 

1212 

62 

3 

HC 

Boulders 

1260 

60 

2 

HC 

Cobbles  (Stones)  -  Shell 

3261 

1261 

61 

3 

HC 

Pebbles  -  Shell 

3262 

1262 

62 

3 

HC 

Granules 

3263 

1263 

63 

3 

mm 

Sand  -  Silt  -  Clay  -  Shell 

3264 

1212 

62 

16 

mum 

Gravel  -  Sand  -  Shell 

3265 

1261 

61 

3 

HC 

Shell 

3266 

1262 

62 

3 

HC 

Rock  -  Gravel  -  Sand  -  Shell 

3267 

1261 

61 

3 

HC 

Sand  -  Shell 

3268 

1213 

63 

4 

HC 

Rock  -  Gravel 

3270 

1201 

1 

2 

HC 

Rock  -  Coral 

3271 

1201 

1 

1 

HC 

Rock  -  Sand 

3272 

1201 

1 

2 

HC 

Rock  -  Mud 

3273 

1 

2 

HC 

Mud  -  Shell 

3274 

1214 

63 

10 

r  HC 

Gravel  -  Sand 

3275 

1213 

63 

4 

HC 

Gravel  -  Mud 

3276 

1214 

63 

10 

HC 

Clayey  Sand  -  Shell 

3277 

1220 

21 

7 

HC 

Soft  Mud 

3278 

1206 

50 

21 

HC 

Hard  Mud 

3279 

1206 

50 

19 

HC 

Silty  Sand  -  Shell 

3281 

1221 

21 

5 
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Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

HC 

Gravelly  Sand  -  Shell 

3282 

1213 

63 

4 

HC 

Medium  Sand  -  Shell 

3283 

1220 

21 

5 

HC 

Fine  Sand  -  Shell 

3284 

1220 

21 

7 

HC 

Sandy  Gravel  -  Shell 

3285 

1213 

63 

4 

HC 

Clayey  Silt  -  Shell 

3286 

1230 

21 

16 

HC 

Silt  -  Shell 

3287 

21 

16 

HC 

Silty  Gravel  -  Shell 

3288 

1214 

63 

10 

HC 

Sandy  Silt  -  Shell 

3289 

1213 

21 

8 

HC 

Sand  -  Mud 

3294 

1209 

24 

17 

HC 

Muddy  Sand 

3295 

1209 

24 

12 

HC 

Sandy  Mud 

3296 

1204 

25 

18 

HC 

Mud 

3297 

1205 

52 

18 

HC 

Clay  -  Shell 

3298 

1240 

21 

10 

HC 

Stiff  Mud 

3299 

1206 

50 

19 

HT 

Rock 

3301 

1101 

1 

2 

HT 

Sand 

3302 

1102 

23 

9 

HT 

Silty  Sand 

3303 

1103 

25 

11 

HT 

Sandy  Silt 

3304 

1104 

25 

16 

HT 

Silt 

3305 

1105 

52 

18 

HT 

Clayey  Silt 

3306 

1106 

50 

19 

HT 

Silty  Clay 

3307 

1107 

50 

22 

HT 

Clay 

3308 

1108 

8 

23 

HT 

Sand  -  Silt  -  Clay 

3309 

1109 

24 

17 

HT 

Gravel 

3312 

1112 

62 

3 

HT 

Sandy  Gravel 

3313 

1113 

63 

4 

HT 

Silty  Gravel 

3314 

1114 

63 

10 

HT 

Muddy  Sandy  Gravel 

3315 

1114 

63 

6 

HT 

Clayey  Gravel 

3316 

1116 

63 

10 

HT 

Muddy  Gravel 

3317 

1116 

63 

10 

HT 

Gravelly  Muddy  Sand 

3318 

1120 

63 

8 

HT 

Gravel  -  Silty  Sand 

3319 

1113 

63 

8 

HT 

Gravelly  Sand 

3320 

1120 

21 

7 

HT 

Very  Coarse  Sand 

3321 

1121 

21 

5 

HT 

Coarse  Sand 

3322 

1122 

22 

7 

HT 

Medium  Sand 

3323 

1123 

23 

9 

HT 

Fine  Sand 

3324 

1124 

24 

11 

HT 

Very  Fine  Sand 

3325 

1125 

25 

13 

HT 

Clayey  Sand 

3326 

1126 

25 

14 

HT 

Gravel  -  Shell 

3328 

1112 

62 

3 

HT 

Gravelly  Silt 

3330 

1130 

21 

16 

HT 

Gravelly  Silt  -  Shell 

3331 

1130 

21 

16 

HT 

Gravelly  Sandy  Silt 

3332 

1130 

21 

16 

HT 

Gravelly  Mud 

3333 

1130 

21 

16 

HT 

Gravel  -  Sand  -  Mud 

3334 

1113 

63 

6 

HT 

Rock  -  Sand  -  Mud 

3335 

1161 

61 

3 

HT 

Rock  -  Gravel  -  Mud 

3336 

1161 

61 

3 

HT 

Rock  -  Gravel  -  Sand 

3337 

1161 

61 

3 
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Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

HT 

Rock  -  Gravel  -  Sand  -  Mud 

3338 

1161 

61 

3 

HT 

Gravelly  Clay 

3340 

1140 

21 

21 

HT 

Sand  -  Clay  -  Shell 

3341 

1102 

23 

9 

HT 

Sandy  Clay 

3342 

1142 

25 

20 

HT 

Very  Fine  Silt 

3350 

1150 

50 

21 

HT 

Fine  Silt 

3351 

1151 

51 

19 

HT 

Medium  Silt 

3352 

1152 

52 

17 

HT 

Coarse  Silt 

3353 

1153 

53 

15 

HT 

Rough  rock 

3354 

1101 

1 

1 

HT 

Mud  over  Rock 

3355 

1101 

1 

2 

HT 

Silty  Clay  -  Shell 

3356 

1130 

21 

16 

HT 

Boulders 

3360 

1160 

60 

2 

HT 

3361 

1161 

61 

3 

HT 

Pebbles  -  Shell 

3362 

1162 

62 

3 

E3 

Granules 

3363 

1163 

63 

3 

Sand  -  Silt  -  Clay  -  Shell 

3364 

1112 

62 

16 

HT 

Gravel  -  Sand  -  Shell 

3365 

1161 

61 

3 

HT 

Shell 

3366 

1162 

62 

3 

Rock  -  Gravel  -  Sand  -  Shell 

3367 

1161 

61 

3 

Sand  -  Shell 

3368 

1113 

63 

4 

HT 

Rock  -  Gravel 

3370 

1101 

1 

2 

HT 

Rock  -  Coral 

3371 

1101 

1 

1 

HT 

Rock  -  Sand 

3372 

1101 

1 

2 

HT 

Rock  -  Mud 

3373 

1101 

1 

2 

HT 

Mud  -  Shell 

3374 

1114 

63 

10 

HT 

Gravel  -  Sand 

3375 

1113 

63 

4 

HT 

Mud  -  Gravel 

3376 

1114 

63 

10 

HT 

Clayey  Sand  -  Shell 

3377 

1120 

21 

7 

HT 

Soft  Mud 

3378 

1106 

50 

21 

HT 

Hard  Mud 

3379 

1106 

50 

19 

HT 

Silty  Sand  -  Shell 

3381 

1121 

21 

5 

HT 

Gravelly  Sand  -  Shell 

3382 

1113 

63 

4 

Medium  Sand  -  Shell 

3383 

1120 

21 

5 

Fine  Sand  -  Shell 

3384 

1120 

21 

7 

HT 

Sandy  Gravel  -  Shell 

3385 

1113 

63 

4 

HT 

Clayey  Silt  -  Shell 

3386 

1130 

21 

16 

HT 

Silt  -  Shell 

3387 

1130 

21 

16 

HT 

Silty  Gravel  -  Shell 

3388 

1114 

63 

HT 

Sandy  Silt  -  Shell 

3389 

1113 

21 

8 

HT 

Sand  -  Mud 

3394 

1109 

24 

17 

HT 

Muddy  Sand 

3395 

1109 

24 

12 

Sandy  Mud 

3396 

1104 

25 

18 

HT 

Mud 

3397 

1105 

52 

18 

HT 

Stiff  Mud 

3399 

1106 

50 

19 

HV 

Clayey  Gravel 

3816 

1116 

63 

17 

HV 

Rock  -  Gravel  -  Sand  -  Shell 

3867 

1161 

61 

3 

HV 

Muddy  Sand 

3895 

1109 

24 

12 
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Table  22  (cont’d.).  Sediment  Type  Definitions  from  the  Enhanced  Database. 

Source:  NAVO  (2006). 


Type 

Enhanced  Categories 
Description 

Enhanced 

Standard 

Reduced 

HFEVA 

HV 

Sandy  Mud 

3896 

1104 

25 

18 

HV 

Mud 

3897 

1105 

52 

18 

HS 

Ooze 

3510 

1108 

8 

23 

N 

Oceanic  Rock  Outcrops 

4000 

1101 

1 

1 

N 

Continental  Rock  Outcrops 

5000 

1101 

1 

2 

N 

Hard  Bottom 

5069 

1101 

1 

3 

A 

Marsh 

6845 

1109 

24 

18 

A 

Mangrove 

6850 

1109 

24 

11 

A 

Intratidal 

6855 

1109 

24 

9 

A 

Supratidal  Zone 

6860 

6860 

21 

999 

A 

Salt  Dome 

6865 

6865 

61 

3 

A 

Gypsum 

6870 

6870 

61 

3 

A 

Peat 

6875 

1102 

23 

12 

A 

Rock  Outcrop 

6880 

1101 

1 

2 

A 

Rubble 

6885 

1162 

61 

1 

A 

Manmade  Features 

6890 

1101 

1 

i 

A  Anomalous 

C  Calcareous 

H  Hemipelagic 

N  Non-Depositional 

L  Land 

O  No  Data 

P  Pelagic 

S  Siliceous 

T  Terrigenous 

U  Turbiditic 

V  Volcanic 
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APPENDIX  B.  SPICINESS  CALCULATION  FUNCTION 


The  following  equation  and  coefficients  for  calculating  spiciness  (Table  Bl) 
are  from  Flament  (1986). 

5)' 

i  =0  i—O 


Table  Bl:  Coefficients,  bn. 

of  the  polynomial  expression  for  spiciness. 

0 

1 

j 

2 

3 

4 

0 

”6 

7.7442e-l 

-5.85e-3 

-9.84e-4 

-2.06e-4 

1 

5.1655e-2 

2.034e-3 

-2.742e-4 

-8.5e-6 

1.36e-5 

i  2 

6.64783e-3 

-2.468 le-4 

-1.428c-5 

3.337e-5 

7.894e-6 

3 

-5.4023e-5 

7.326e-6 

7.0036e-6 

-3.041 2e-6 

-1.0853e-6 

4 

3.949e-7 

-3.029e-8 

-3.8209e-7 

1.0012e-7 

4.7133e-8 

5 

-6.36e-10 

-1.309c-9 

6.048e-9 

-1.1409e-9 

-6.676e-10 

A  Matlab  version  of  Flament’s  algorithm  is  listed  below  as  SPICE.M: 

SPICE.M:  Spiciness  Calculation  Function 

function  Tau  =  spice(S, Theta) 

%  SPICE  -  'Spiciness1,  an  oceanographic  variable  for 
%  characterization  of  intrusions  and  water  masses. 

% 

%  Use  as:  Tau  =  spice(S.Tbeta) 

%  Inputs:  S  =  Salinity  (psu) 

%  Theta  =  Potential  temperature  (°C) 
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%  Output:  Tau  =  spiciness 
% 

%  Brian  Schlining;  07  Oct  1997 
%  Algorithm  from  Flamcnt,  1985. 

%  Table  of  Coefficients  b(i,j)  of  polynomial  expression 
b  =  [  0  7.7442e-l  -5.85e-3  -9.84e-4  -2.06e-4;... 

5.1655e-2  2.034e-3  -2.742e-4  -8.5e-6  1.36e-5;... 
6.64783e-3  -2.468 le-4  -1.428e-5  3.337e-5  7.894e-6;... 
-5.4023e-5  7.326e-6  7.0036e-6-3.0412e-6-1.0853e-6;... 
3.949e-7  -3.029e-8  -3.8209e-7  1.0012e-7  4.7133e-8;... 
-6.36e-10  -1.309e-9  6.048e-9  -1.1409e-9-6.676e-10]; 

%  Polynomial  calculation  of  spiciness 
Tau  =  b(l,l).*Theta.A0.*(S-35).A0  + ... 
b(l  ,2).  *Theta.A0.*(S-35). A 1  +  ... 
b(l,3).*Theta.A0.*(S-35).A2  + ... 
b(l,4).*Theta.A0.*(S-35).A3  + ... 
b(l,5).*Theta.A0.*(S-35).A4  +  ... 
b(2,l).*Theta.Al.*(S-35).A0  + ... 
b(2,2).*Theta.Al.*(S-35).Al  +  ... 
b(2,3).*Theta  A1.*(S-35).A2  + ... 
b(2,4).*Theta.Al.*(S-35).A3  +  ... 
b(2,5).*Theta  Al.*(S-35)  A4  + ... 
b(3,l).*Theta.A2.*(S-35).A0  +  ... 
b(3,2).*Theta  A2.*(S-35).A1  + ... 
b(3,3).*Theta  A2.*(S-35).A2  + ... 
b(3.4).*Theta  A2.*(S-35).A3  + ... 
b(3,5).*Theta  A2.*(S-35)  A4  + ... 
b(4,l).*Theta  A3.*(S-35).A0  + ... 
b(4,2).*Theta.A3.*(S-35).Al  + ... 
b(4,3).*Theta.A3.*(S-35)  A2  + ... 
b  (4,4).  *Theta.  A3.*(S-35).A3  + ... 
b(4,5).*Theta  A3.*(S-35).M  + ... 
b(5,l).*Theta  A4.*(S-35)  A0  + ... 
b(5,2).*Theta  A4.*(S-35).A1  + ... 
b(53).*Theta.A4.*(S-35).A2  +  ... 
b(5,4).  ^Theta.  A4.  *(S-35).  A3  + ... 
b(5,5).*Theta.A4.*($-35).A4  + ... 
b(6,l).*Theta.A5.*(S-35)  A0  + ... 
b(6,2).*Theta.A5.*(S-35)  A1  + ... 
b(6,3).*Theta  A5.*(S-35).A2  + ... 
b(6,4).*Theta  A5.*(S-35)  A3  + ... 
b(6,5).  *Theta.A5 .  *(S-35).  A4: 
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Latitude  Latitude  Latitude  Latitude  Latitude  Latitude 


APPENDIX  C.  ADDITIONAL  GDEM  MONTHLY  SPICE  FIGURES 


GDEM  Monthly  Spice  at  400  m 
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Latitude  Latitude  Latitude  Latitude  Latitude  Latitude 


GDEM  Monthly  Spice  at  1000  m 
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APPENDIX  D.  ADDITIONAL  SMG-WOD  MONTHLY  MEAN 

SOUND  SPEED  FIGURES 


SMG-WOD  Monthly  Mean  Sound  Speed  at  400  m 
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SMG-WOD  Monthly  Mean  Sound  Speed  at  1000  m 
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Latitude  Latitude  Latitude  Latitude  Latitude  Latitude 


APPENDIX  E.  ADDITIONAL  GDEM  MONTHLY  SOUND  SPEED 

FIGURES 


GDEM  Monthly  Sound  Speed  at  400  m 
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Latitude  Latitude  Latitude  Latitude  Latitude  Latitude 


GDEM  Monthly  Sound  Speed  at  1000  m 
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